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Abstract 

This paper gives detailed information on the structure of the infrared singularities for the 
process e^e^ 3 jets at next-to-next-to-leading order in perturbation theory. Particular em- 
phasis is put on singularities associated to soft gluons. The knowledge of the singularity 
structure allows the construction of appropriate subtraction terms, which in turn can be im- 
plemented into a numerical Monte Carlo program. 



1 Introduction 



The process e^e^ 3 jets is of particular interest for the extraction of the strong coupling a^. 
Observables related to three-jet events are well suited for this task because the leading term 
in a perturbative calculation of three-jet observables is already proportional to the strong cou- 
pling. An accurate and precise extraction of tty from the experimental measured observables 
requires also a precise theoretical prediction for these observables. Infrared-safe observables 
can be computed in perturbation theory and higher precision is obtained by including the next- 
higher term in the perturbative series. Infrared-safe observables related to e+e — > 3 jets have in 
the past been calculated at leading-order (LO) and next-to-leading order (NLO) [1-6] in QCD. 
The NLO electro- weak corrections have been calculated recently in [7]. This paper reports on 
the next-to-next-to-leading order (NNLO) QCD calculation for e+e^ 3 jets. This is a highly 
non-trivial calculation which triggered research in many directions. The relevant tree-, one- and 
two-loop-amplitudes have been calculated in refs. [5,8-19]. For the computation of the two-loop 
amplitudes new integration techniques have been developed [20-23] together with numerical 
routines for the evaluation of poly logarithms [24-26], which occur in the integral functions. The 
inclusion of the third term of the perturbative expansion into a numerical program, which allows 
the calculation of an arbitrary infrared-safe three-jet observable is very challenging. Most of the 
complications are related to infrared singularities, which occur in intermediate steps. Methods to 
handle these singularities at NNLO have been discussed in refs. [27-44], generalising ideas from 
NLO [5,45-51]. Unlike other processes like e+e" — > 2 jets, pp ^ H or Drell-Yan, which have 
been calculated previously to next-to-next-to-leading order [52-68] and which involve only two 
hard coloured partons, the process e^e~ —>■ 3 jets is the first process calculated at NNLO with 
three hard coloured partons. 

Here new complications related to soft gluons occur. This paper documents the singularity 
structure of — > 3 jets and gives a detailed discussion of the singularities related to soft 
gluons. Another group published results on the infrared structure of e~^e~ 3 jets at NNLO 
earlier on in [69], but omitted certain subtraction terms related to soft gluons. In the calculation 
presented here the methods used are in many parts similar to the ones used in [69] and it should 
be mentioned that the authors of [69] made major contributions to the development of these 
methods [33-37]. The purpose of this paper is to present a complete set of subtraction terms for 
the process e+e 3 jets. Meanwhile the other group has added the missing subtraction terms 
in the colour factors A'^^ and A'^. 

In general the subtraction terms are not unique, since arbitrary finite terms can be added 
to them. There are a few minor points, where the subtraction terms of this paper differ from 
the corrected subtraction terms of ref. [69]. First of all, ref. [69] introduces "dipole" subtraction 
terms by dividing a three-parton antenna into two dipoles. The subtraction terms in this paper are 
based entirely on complete antennas. Secondly, ref. [69] symmetrises in the quark momenta, thus 
eliminating contributions anti- symmetric in the quark momenta. In this paper I give subtraction 
terms which approximate the double-real emission contribution in all singular limits without the 
need for symmetrisation in the quark momenta. Thirdly, ref. [69] uses in the colour factor N/Nc in 
the NLO subtraction terms for a coUinear quark pair always an antenna with a gluon as spectator. 
As a consequence no soft subtraction term is needed with this choice. In this paper the NLO 
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subtraction terms are the canonical ones, obtained from the colour flow of the amplitudes. In this 
case a soft subtraction term occurs also in the colour factor N/Nc. In general the occurrence of 
soft subtraction terms is related to the specific choice of the NLO subtraction terms. In this work 
the subtraction terms are developed from the starting point, that the NLO subtraction terms are 
given by the canonical choice. A second possibility, which completely avoids soft subtraction 
terms, is to modify the NLO subtraction terms appropriately. This possibility has been discussed 
recently in ref. [70]. 

The subtraction terms used in this work are based on spin-averaged antenna functions. There- 
fore they do not capture phase-space singularities associated to spin-correlations and are there- 
fore not local subtraction terms. These remaining singularities are rather mild and treated by a 
slicing procedure. Therefore the complete treatment of all infrared singularities for the process 
e~^e~ 3 jets is based on a hybrid method of subtraction and slicing. 

This paper is necessarily rather technical. It was written for the following reason: Under- 
standing the infrared singularity structure of the process e^e — > 3 jets at NNLO is not only 
relevant to three-jet observables in electron-positron annihilation, but to all processes with three 
or more hard partons, if calculated at NNLO. The knowledge of the singularity structure allows 
the construction of appropriate subtraction terms and I give a full account of the subtraction 
terms used in the numerical program "Mercutio2" [71,72]. Numerical results for jet rates and 
event shapes related to the process e^e~ — > 3 jets have been presented in [69,71,73-75]. The 
numerical results from the program "Mercutio2" are presented in a companion paper [76]. 

This paper is organised as follows: Section 2 introduces the notation and recalls the most 
important facts of the relevant amplitudes. Section 3 discusses the general features of the method 
used to handle the infrared divergences. The explicit subtraction terms are given in section 4. 
Section 5 shows explicitly the cancellation of all singularities. Section 6 contains a summary 
and the conclusions. In appendix A an alternative form for the soft subtraction term is given. 
Appendix B lists the explicit forms of the various terms in the singular limits. 



2 Notation and conventions 



2.1 Cross sections and amplitudes 

The master formula to calculate an observable at an collider with no initial-state hadrons (e.g. an 
electron-positron collider) is given by 

11 1 f ' 

^^^'^^ = (OT ^UfOT ^^"^"'^ ^^1' -.Pn^quqi) L \M\ (1) 

2K[s) (2/i + 1 ) (2/2 + 1 ) T heUcity 

where q\ and q2 are the momenta of the initial-state particles, 2K{s) — 2s is the flux factor and 
s— (^1 -1-^2)^ is the center-of-mass energy squared. Kinematical invariants will be denoted by 

2 2 

Sij={pi + Pj) , Sijk= (pi + Pj+Pk) , etc. (2) 

The factors 1 / (27i + 1) and 1 / (2J2 + 1) correspond to an averaging over the initial helicities. 
d^n is the invariant phase space measure for n final state particles and On ^ {pi,...,p„,qi,q2) is 
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the observable, evaluated with a configuration depending on n final state partons and two initial 
state particles. The index j indicates that the leading order contribution depends on j final-state 
partons. The observable has to be infrared safe, in particular this implies that in single and double 
unresolved limits we must have 



oiii{Ph-,Pn+l,qhq2) 
0l/^2iPh-,Pn+2,qhq2) 



On\p\,...,Pf^,q\,q2) for single unresolved limits, 

{p'\i---i p'nidx 1 4i) for double unresolved limits. (3) 



J?„ is the amplitude with n final-state partons. At NNLO we need the following perturbative 
expansions of the amplitudes: 

1^/ = (^f + ^W^^i*^)) + (^f + ^^^^ 

2_ ^(0) *^(0) ^^,(0) 



|-^n+l| — ^n-ll + i "^n+l ^n+\ + "^n+l "^n+l) ' 



^(0) *_^(0) 



(4) 



Here sin denotes an amplitude with n final-state partons and / loops. We rewrite the master 

(5) 



formula eq. (1) symbolically as 

n •' 

and the LO, NLO and NNLO contribution as 

I^qU)\jNlo 

I^O^j))NNLO 



On do II , 



(6) 



-> 3 jets the colour decomposition reads 



The LO, NLO and NNLO contributions may be decomposed into gauge-invariant colour pieces. 
For the process e^e~ 

{^0^^))NNLO 



2 
1 



_i_ AT- Al lr> (3) \NNLO ,^f/r, (3) \NNLO , / ^ (3) \NNLO 
+NfNc{0^ '}nf +j;^{0^ ')nfsc +^f{0^')nfnf 



(7) 



where Nc denotes the number of colours and Nf the number of light quark flavours. I use the 
following normalisation for the colour matrices: 



Tr T^'T^ — —5"^ 
2 



(8) 
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Table 1: Contributions at NNLO of the various amplitudes to the individual colour factors. 



The factors 1/2, 1/4 and 1 /8 in eq. (7) are a consequence of this normalisation. In addition to the 
terms shown in eq. (7), there are in the colour decomposition singlet contributions, which arise 
from interference terms of amplitudes, where the electro-weak boson couples to two different 
fermion lines. These singlet contributions are expected to be numerically small [18,77,78] and 
neglected in the present calculation. 

The computation of the NNLO correction for the process e+e" 3 jets requires the knowl- 
edge of the amplitudes for the three-parton final state e'^e~ — > qqg up to two-loops [18, 19], 
the amplitudes of the four-parton final states e'^e~ — > qqgg and e^e~ — > qqq'q' up to one-loop 
[13-16] and the five-parton final states ^Iggg and e~^e~ qqq'q' g at tree level [8,9,79]. 

The one-loop amplitude e+e^ ggg enters only the NNLO term as a quark singlet contribution 
and is therefore neglected. I will label the momenta of the particles of the amplitudes as follows: 

-^4!L^(0'1'2,3) = Jl[^l^^^{qo,gi,g2,q3,e4,e+), 

^4%'q'q(^^'^^^^^) = ^i%'^'qi^o,q[,q'2,q3,e4,e+), 

■^5!L^9(0'1'2,3,4) = Sl^5lggg^{qo,ghg2,g3,q4,e^,e+), 

-^5%W.(0'1'2'3,4) = 4%'gq'q(^0,q'ug2,q'3,q4,e;,et). (9) 

The ordering of the external particles follows the conventions of [80]. I use the convention that all 
particle momenta are outgoing. Therefore in eq. (9) corresponds to an outgoing electron with 
momentum pi, or equivalently to an incoming positron with momentum —pi. Not all amplitudes 
contribute to all colour factors. Table 1 shows the contributions at NNLO to the individual colour 
factors. 



2.2 Colour decomposition of the amplitudes 

For the discussion of intermediate results it is useful to have the colour decomposition of all the 
relevant amplitudes. We start with the amplitudes for the three-parton final state e~^e~ qqg. 
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The colour structure of the tree amplitude, the one-loop amplitude and the two-loop amplitude is 



^3^\^0,ghq2) = e^gco{T^)^QAf\qo,gi,q2), 

4^\qo,8uq2) = e^8Co{T')^^{^y4\qo,guq2). (10) 

Here, e denotes the electro-magnetic coupling, g the strong coupling and as usual we set a = 
e^/ (4%) and = ^2/(471). cq denotes a factor from the electro-magnetic coupling: 

co = -!2" + vVfPz(*), 2'z(^) = • (11) 

s — + li. zmz 

The electron - positron pair can either annihilate into a photon or a Z-boson. The first term in 
the expression for cq corresponds to an intermediate photon, whereas the last term corresponds 
to a Z-boson. The left- and right-handed couplings of the Z-boson to fermions are given by 

J l(-Qfsm^Bw f -QfsinQw .... 
smQwcosQw cosQw 

where and are the charge and the third component of the weak isospin of the fermion. If 

the anti-quark qo in eq. (10) has colour index Jq, the quark q2 has colour index (2 and the gluon 
gi has colour index ai, we abbreviate the notation for the colour matrix by 

The one-loop amplitude Jl^^^ has the colour decomposition 

^(1) - +1^(1) (14) 

-^slL- -^3 ax represent contributions, where the electro-weak gauge boson couples to a closed 
fermion loop with a vector or axial-vector coupling, respectively. Due to Furry's theorem (i.e. 

charge conjugation) the amplitudes -^^3^!^^ vanish. Furthermore, since the u, d, s, c quark masses 
may be neglected only the t, b quark pair survives an isodoublet cancellation in the fermion loop 
and contributes to ^^Ij^- Explicit results for the amplitude J^g^^ can be found in refs. [5, 14]. The 

' (2) 

two-loop amplitude Sl^ has the colour decomposition 

^3 - 4 -^3,/c^4'^3,.-c+4^2'^3,,s-.vc+ 4 '^3,n/^ 4^^'^3,n/ic+ 4 '^3,n/n/ 

+ ^^(2) +^^(2) , ^ (2) ,^.(2) 
' 4 3,vec,lc ' 4 3,ax,lc ' 4^ 3,vec,sc ' 4^ 3,ax,sc' 
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(2) (2) (2) (2) 

■^3,vec,/c' ■^3,vec,.c' Hax,ic ^3,ax,sc represent again contributions, where the electro-weak 
gauge boson couples to a closed fermion loop with a vector or axial-vector coupling. These 
contributions are expected to be small and neglected in the present analysis. Explicit results for 

(2) 

the amplitude Jl^^ can be found in refs. [18, 19]. 

We now turn to the colour decomposition of the amplitudes of the four-parton final states. 
We first consider the decomposition of the amplitudes for the sub-process e'^e~ — > qqgg. The 
colour decomposition of the Bom amplitude is 

4''\^0,ghg2,q3)=e^g^co[{T^T')^QAf\qo,gi,g2,q3) + {T'T^)^QAf\q^^ 
The one-loop amplitude is first decomposed into partial amplitudes 
K W0,ghg2,q3) =e g Co 

A^c((r2rl)3oAy(^0,gl,g2,?3) + (rlr2)3oAi|;(<7o,^2,gl,?3)) 
+ l^'%oA'i^l{qo,ghg2,q3)+(^{T^T')3o + {T'T^^^^ 

+ {(T'T'ho-^f%oyiijqo,gl,g2,q3) 



+ {('^'T'ho-:^f%oyi\^qo,g2,gu^^^ 



.(17) 



In this formula the one-loop leading-colour partial amplitude is denoted by A^^l and A^^l denotes 

a sub-leading-colour contribution. The partial amplitudes ^4']^^,, A^^] ^ and a''^^ ^ represent the 
contributions from a photon or Z coupling to a fermion loop through a vector or axial-vector cou- 
pling. These contributions are small and neglected in the further analysis. The partial amplitude 
A^^j may be further decomposed: 

^i^J(^0,gl,g2,?3) = 

^i^l\zc(^0,gl,g2,?3) + ^Ay,„/(?0,gl,g2,?3) + -^Ai|{^^^(^0,gl,g2,?3). (18) 

Explicit results for the amplitude ^4^\qo,gi,g2,q3) can be found in ref. [14]. 

We now consider the colour decomposition of the amplitudes for the four-parton final state 
qqq'q'. The colour decomposition of the tree amplitude reads 

^i^\q0,qi,q2,q3) = I ^ (^S3i520- -^530521^ (^0,^1,^2: ?3) 

-5//av^ (^830821 - -^531520^ Xi*^^ (^1,^0, ^2, ^3) 
X4°^(^0,^i,?2>?3) = Co(3)Af^(^o,^i,?2>?3)+Co(2)A^°(^l,^0,?3,?2)- (19) 
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Here co{q) denotes the electro- weak coupling for quark q. The decomposition of the one-loop 
amplitude is 

R^^\qo,q\,q2,q3) = \ ^0 , ^ i , ^2 ^ ^3 ) - 5//avXi^ ^ (<? 1 , '^O , > ^3 ) , 

X4^H^0,^i,?2>?3) = eV (^) 2 (^c53lS20X4!l(^0,^i,?2>?3) +53052ixi2(^0,^i,?2>?3) 

14} (^0,^1 , ?2> ?3) = CO (3)A^^J (^0, q'l , ?2> ?3) + co(2)Ai|} (^1 , ^0, ?3 , ^2) > 

xi|2(^o,^i,^2,^3) = co(3)Ai|](^o,^i,^25^3)+co(2)Ai|^(^i,,?o,^3,^2)- (20) 

The function %4^]^ arises from a fermion triangle graph with an axial coupling of the Z and is 

neglected in the further analysis. The partial amplitudes A4^j and A4^2 can be further decomposed 
according to 

^4!l(^0,^l,?2>?3) = 

4!l,/c(^0,^i,?2,?3) + ^^A^4lnA^0,q'i,q2,q3) + ^A[]l^^{qo,q[,q'2,q3), 

4!2,/c(^o,^i,?2,?3) + jf^A^4,2,nf{qo,q'uq2,q3) + ■^^4lsMo,q'uq2,q3)- (21) 

Explicit results for the amplitude R^^^ {qo, q[ , ?3) found in ref. [13]. 

Finally let us discuss the colour decomposition of the amplitudes related to the five-parton 
final state. The colour decomposition of the tree-level amplitude for e~^e~ qqggg is 

^Piqo,8hg2,g3,q4) = 

{T^T^T')4oAf\qo,gi,g2,g3,q4) + {T^T^T')4oAf\qo,gi,g3,g2,q4) 
+ {T'T^T^)4QAf\qo,g2,g3,gl,q4) + {T'T'T^)4oAf\qQ,g2,gi,g3,q4) 
HT^T'T^)4oAf\qo,g3,gi,g2,q4) + iT'T^T^)4oAf\qo,g3,g2,gl,q4). (22) 

Explicit results for the amplitude ^5*^^ (^0,^1, ^2, ^3, ^4) can be found in ref. [8,79]. 
The colour structure of the tree-level amplitude for — > qqq'q'g is given by 

■^f^(^0,^i,g2,?3>?4) = ]^^A0T^\Dx{qQ,q\,g2,q'3,q4) + ]^h\T^D2{qQ,q\,g2,q'3^q4) (23) 

- o S41 7^30^3 (^0, ^1 , .?2, ?3 ' q^) - o ^30^41^4 (^0, qi,g2,q3, q^) ■ 
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The functions Dj{qo, q\^g2i ^3,^4) can be further decomposed according to 

^1(^0,^1, ^2, ?3:?4) = Co(4)5i(<?o,<?i,^2,^35'?4) +Co(3)52(^i,^0,.?2,^4,'?3) (24) 

+^flavTT (<^0(4)53(^l,^0,g2,?3'?4) + C0(3)54(^0,^l,g2,?4,?3)) ' 

Diiqa, q'l , g2 , ?3 > ?4 ) = CO (4)52 (^0 , ^ 1 , g2 , ?3 > ?4 ) + CO ( 3 )5 1 1 , ^0 , ,?2 , ?4 , ?3 ) 

+^flavj^ {co{4-)B4{q\,qo,g2,q3,q4) + Co{3)B3{qo,q[,g2,q4,q3)) , 

^3 {qo , q'l , g2 , ?3 , ?4) = CO (4) ^83 {qo , ^'j , ^2 , ^3 , ?4) + Co (3 ) (<?i , ^0 , g2 , ?4 , ?3 ) 

+5/Zav (co(4)fil(^i,^0,,?2,<?35'?4)+Co(3)52(^0,^l,^2,?4,'?3)) > 

£>4(^0,^l,g2,?3,?4) = Co(4)^54(^o,^l,g2,?3>?4) + Co(3)^53(^l,^0,g2,?4,?3) 

+5//av (co(4)B2(^l,^0,^2,?3>?4) + Co(3)Bi(^0,?'l,g2,?4,?3)) • 

ExpUcit results for the ampUtude Sl^^^ {qaAi ■> ^2, ^35 ^4) can be found in ref. [8, 79]. 
2.3 Ultraviolet renormalisation of loop amplitudes 

The loop amplitudes have explicit divergences. It is common practise to regulate these diver- 
gences by dimensional regularisation. The origin of these divergences are either related to ultra- 
violet or to infrared singularities. I first discuss ultraviolet divergences. These are removed by re- 
defining the parameters of the theory. For the QCD corrections to the process e^e — > 3 jets with 
massless quarks it is sufficient to renormalise the strong coupling constant. In the MS scheme 
the relation between the bare coupling ao and the renormalised coupling a^(//^) evaluated at the 
renormalisation scale 1^ reads: 



ao = CLsS^^i?^ 



(25) 



where 

5e = (47I)^e-^Y^ (26) 

is the typical phase-space volume factor in D = 4 — 28 dimensions, is Euler's constant, and Po 
and Pi are the first two coefficients of the QCD P-function: 

Po = y Ca - ^ W/, Pi = y cl - ^CATRNf - ACfTrNj, (27) 



with the colour factors 



C,=/V„ C, = ^|^, r„ = i. (28) 
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Let the expansion in the strong coupUng of the unrenonnaUsed ampUtude be 



n-2 



^n,bare = (47iao) ^ 



^(') + 
n,bare ' 



2J \bare + {2%) "'^^''^ ^ ^ 



Then, the renormalised two-loop amplitude can be expressed as 

^(0) /OCA - en /a.9\2 



■^n,ren 1 \ ^-^J ^^^^^ ^ y 271/ V"'* / 



(29) 



(30) 



The relations between the renormalised and the bare amplitudes are given by 



^(0) 

■^n.ren 



n,barei 



- 5-V^^(i) -in-D^k^^^ 



^(2) 
■^n.ren 



-2..4e-(2) „P0e_i.2e-(l) , (" " 2) Po Pl 



2e2 



(0) 

n,bare' 



(31) 



2.4 Infrared structure of loop amplitudes 



Let us now turn to the infrared singularities of loop amplitudes. The infrared pole structure of 
loop amplitudes in the dimensional regularisation parameter £ is well understood [81-85] and 
can be stated explicitly for one- and two-loop amplitudes. If we regard a loop amplitude as a 
vector in colour space, then the infrared poles of the loop amplitude can be expressed through an 
operator acting on this vector in colour space. We start with the one-loop amplitude, which can 
be written as 



(32) 



Here I^^) (e) contains all infrared double and single poles in 1/e and -^^^ is a finite remainder. 
At two-loops, the corresponding formula reads: 



= i(i)(8)^i^)+i(2)(8)^r+^(5,„. 

The one-loop insertion operator I^^) is given by 
where 



and the coefficients T? and y; are 



(33) 



(34) 



(35) 



(36) 
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In general, the colour operators T/Tj give rise to colour correlations. However, the colour struc- 
ture for the Bom amplitude e~^e~ qqg is rather trivial and the colour operators are in this case 
proportional to the identity matrix in colour space: 

Wf = 7-.i-;,f. T,T,<' = T,T,<> = (37) 

The two-loop insertion operator has the form 

,(2)(e) = _l,(i)(e)(i(i)(e) + ^)+,-w.E|z|)(|+^),(i)(2e) + HP), (38) 

where 

^18 6 y 9 

The operator H^^^ is process- and scheme-dependent and for e+e" qqg, it is given by [17, 86] 

1 



«^'^^3'^ = ^mb)^"4K+^^'^+<0^3^ (40) 



where H^^ — H^^ and 

r(2) _ rZr._11.2^409\ ,2 f I. 41 \ / 3. 3 \ I 



<^ = u^-96^+ 864j^^n-i^^-%-T^Jn-2^^+Y-^j^ 

^\4S 216 J N, ^' 



(2) _ fk^^^2^5\.^f_ie_S9\^^^^_Nf 



m^' = + + _ WV^ + N '- + —Ni. (41) 



3 Cancellation of infrared divergences 
3.1 Singular behaviour in phase space 

QCD amplitudes become singular in phase space, when the momenta of one or more external 
particles become degenerate. In perturbative calculations this phenomen occurs first in next- 
to-leading order calculations and is related to the term d<5^^^^. Singularities occur if either two 
particles become coUinear or if one particle becomes soft. The factorisation properties are most 
conveniently discussed through decomposing QCD amplitudes into primitive amplitudes, as it 
was done in section 2.2. Primitive amplitudes are defined by a fixed cyclic ordering of the QCD 
partons, a definite routing of the external fermion lines through the diagram and the particle 
content circulating in the loop. In the soft limit one has the factorisation 

l^^.^nlli-^Pi^Pj^Pk,-) = ^^^f\pi,Pj,PkWn\-,Pi,Pk,-), (42) 

; soft ^ 
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where the eikonal factor is given by 

Eik<»'(p,,P„K) = }iL<EA.}E^, (43) 

The square of the eikonal factor is given by 

J^Eikf (p,-,p;,p,)*Eikf = 4^. (44) 

Soft singularities lead to colour correlations among the primitive amplitudes aI'^^ 

In the coUinear limit one has the factorisation 

= Splitf \pi,pj)Al^\...,p,...). (45) 

where pi and pj are the momenta of two adjacent legs and the sum is over all polarisations. 
Squaring the splitting amplitudes one obtains 

P3%arkM = Splitf (/7,-,/7y)*Splitf (/7,-,/7,-)m(P,V), 

P3%uonM = Y^^'ip^r Splitf\pi,pj)*Sp\itf\pi,pj)E\p,r). (46) 

Xi,Xj 

P^^") (A,, A,') is a tensor in spin space. CoUinear singularities lead therefore to correlations in spin 
space. The spin-averaged splitting functions are obtained by 

i^kquark) = 2 ^ ^3, quark ^) ' 



X 



One-loop amplitudes factorise in single unresolved limits (i.e. soft or coUinear limits) as [28,87- 
92] 

limA^i = i:(sing(°).Ai^) + SingW.Ar). (48) 

In the soft limit, the sum involves only one term, whereas in the coUinear limit the sum is over the 
polarisations of the intermediate parton. As in the previous case, the sum over the polarisations 
leads to correlations in spin space. The one-loop splitting amplitudes have a decomposition into 
primitive splitting amplitudes: 

Sing(i) = Sing(;)^ + ^Singg^ + ^Sing(;i. (49) 
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The one-loop splitting amplitudes have an ultraviolet divergence and are renormalised according 
to 



Singg,„ = 5rV'^Sing(;],,,-^Sing(X.„ (50) 

where is the renormalisation scale. At NNLO we have to consider in addition double un- 
resolved limits related to the term dc^^^2- Generic double unresolved limits are the double soft 
limit, the soft-coUinear limit and the triple coUinear limit. Tree amplitudes factorise in the double 
unresolved limits as [27, 93-98] 

- ESingr-Al^). (51) 

In the double soft limit the sum involves only one term, while in the soft-coUinear limit and in 
the triple coUinear limit the sum is over the polarisations of the intermediate parton. Again, the 
sum over the polarisations leads to correlations in spin space. The double unresolved limit of 
two coUinear pairs factorises in a product of two splitting functions. 

3.2 The subtraction method 

The individual contributions on the r.h.s. of eq. (6) to l^o'^^^)^^^ and (^o^-^^)^^^^ are in general 
infrared divergent, only the sum is finite. However, these contributions live on different phase 
spaces, which prevents a naive Monte Carlo integration approach. To render the individual con- 
tributions finite, one adds and subtracts suitable chosen terms. The NLO contribution is given 
by 

(oO-)r^ = /(o(i\jait-o(^-)orf'^))+|(o(^-)rf + o(^-)orf (52) 

The notation On ^ o da^'^^ is a reminder, that in general the approximation is a sum of terms 

ol^-^orf^) = l^ol^^da^''^ (53) 

and the mapping used to relate the n -I- 1 parton configuration to a n parton configuration differs 
in general for each summand. 

In a similar way, the NNLO contribution is written as 

+ / (oi'-' rfo?' + oi^* o<Ja<«'2' + o"' o<Ja<"') . (54) 
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^cc^+i^ is the NLO subtraction term for (n+ l)-parton configurations, da^'^^ and da^'^^ are 
generic NNLO subtraction terms. It is convenient to rewrite these terms for a process Uke 
e'^e~ — > 3 jets as 

da^n'^^ = daf"*^^' + da'li"""' + da'f - Jaf 

da^l'^'^ = Ja;f''^+t^af''^"'^^-^i<''«''^^-t^<''^'+daf''«''^^. (55) 

Each term has an interpretation. dfX'^'^l'^ is the NLO subtraction term for (n + l)-jet observables, 
containing only three parton tree-level antenna functions. At NNLO there are several new sub- 
traction terms required, each of them with a specific structure. The term da'^""^^'^ contains the 
four-parton tree-level antenna functions. The term da^""^ contains three-parton one-loop antenna 
functions together with tree-level matrix elements and three-parton tree-level antenna functions 
together with one-loop matrix elements. The remaining terms da^^^'''^^^, da^^"^^', doJ"''^^"^* and 
dffof' all contain a product of two three-parton tree-level antenna functions. In Ja'<«''«<e^ and 
^^aimost Qj^g antenna function has five-parton kinematics, while the other antenna has four-parton 
kinematics. The former subtraction term is an approximation to Ja""*''*^, while the latter ap- 
proximates do^^^2 ill almost colour-correlated double unresolved configurations. Almost colour- 
correlated double unresolved configurations correspond to two three-parton tree-level antenna 
functions sharing one common hard radiator. In ddP^"^^ both antennas have four-parton kine- 
matics. The term da^°^' approximates soft gluon singularities. The decomposition of eq. (55) 
applies to the process e~^e~ 3 jets. For processes with more than three hard partons one has 
in addition products of three-parton tree-level antenna functions corresponding to colour discon- 
nected contributions. 



3.3 Antenna functions 

For the subtraction terms we need functions which interpolate between the various singular lim- 
its, such that each limit is exactly approximated once and not over-subtracted. A possibility are 
antenna functions [27, 28, 33-38, 53, 99]. I will largely follow the notation of ref. [37]. In a 
few smaller points the notation differs and this paragraph defines the notation used in this article 
for the antenna functions. The small modifications simplify the presentation of the subtraction 
terms. Generically the unintegrated antenna functions are denoted by 

X§ii,j,k), Xiii,j,k), X^ii,j,k,l), (56) 

where denotes a three-parton tree-level antenna function, denotes a three-parton one- 
loop antenna function and X^ denotes a four-parton tree-level antenna function. A pictorial 
representation for these basic antenna functions is shown in fig. 1. At NNLO we also have 
iterations of the three-parton tree-level antenna functions. For three hard coloured partons the 
two generic cases shown in fig. 2 occur. They are labelled 'colour-connected' and 'almost colour- 
connected' . For four or more hard partons also colour disconnected configurations occur. These 
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single unresolved 



one-loop unresolved 



double unresolved 



Figure 1 : The basic antenna functions occurring in an NNLO calculation: The three parton tree- 
level antenna function (left), the three parton one-loop antenna function (middle) and the 
four-parton tree-level antenna functions X^ (right). 




colour connected 




almost colour connected 



Figure 2: At NNLO also the iteration of the three parton tree-level antenna functions X9 occur. 



are rather trivial and not shown in fig. 2. 

Letters from the beginning of the alphabet are used to denote specific antenna functions: 



The relevant four-parton antenna functions are 

^4('V' jg^ h 

^^{iq, jq, kq, Iq 
C4 {iq, jq-, kq, Iq 



gj+gk 
gj+gk 



Gi+gj + Gk, 
Qi+gj + gk, 
qi + qj + qk, 
gi+gj + gk, 
gi + qj + qk- 



qj + qj qi + gj+gk + qi, 
qj + qf qi + q'j + q'k + 
qj + qf^ qi + qj + qk + qi, 



(57) 
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^4 (iqi Jg ikgjg): qj+gi^ qi + gj +gk+gh 

E^iiq, jq, kq, Ig) I ^- + gj ^ qi + qj + qk + gl- (58) 

A subscript "sc" or "nf" is used to indicate contributions sub-leading in colour or proportional 
to the number of light flavours Nf. (In ref. [37] a tilde and a hat are used to denote these contri- 
butions.) The leading colour one-loop three-parton antenna function is denoted with a subscript 
"Ic". For the four-parton tree-level antenna functions no extra subscript is used for the contribu- 
tion leading in colour. For the three-parton one-loop antenna functions we therefore have 

^3,lci^qJgJ^q)j ^\,nf{^qijg-'^q)i ^\,sc{^qijgi^q)i 
^\,lc{^qi jgi^g) ' ^\,nf{hi jgi ^g)i 
^3,lci^^jq^^q)^ ^3,nfi^q^jq^^q)^ ^3,sci^T jqT^q)^ 
^3,lcih^ jg^ ^g)^ ^3,nf{h^ ^g)^ 

^3,tcih^Jq^^q)^ ^3,nfih^Jq^^q)^ ^3,sAh^jq^^q)- (^^) 
We denote by the combination 

The minus sign in the sub-leading colour term is convention. For the four-parton tree-level 
antenna functions we have 

^4('(?5 Jg7 kg, Iq) , A^ .^^{iq, jg,kg, Iq) , 

E4{iq,jq,kq,lg), E^^^^[iq, jq,kg, Iq) . (61) 

The arguments of the antenna function ^^{iq,jq,kg, Iq) are labelled according to the singularity 
structure. The four-parton tree-level antenna functions B'^, and D4 have only one colour 
contribution. 1 will assume that the antenna functions are already ordered. Therefore the antenna 
functions contain only singularities corresponding to the order of the arguments. In particular the 
first and the last argument of the ordered antenna functions indicate the hard partons. This is the 
major difference with the notation of ref. [37]. This concerns the three-parton antenna functions 

D%,jg,kg), Fl{igjg,kg), / = 1,2 (62) 

and the four-parton antenna functions 

'^4,jc('?'-/g'^g'^9)' ^^i}qT jg^kgilg), E^[iq, jq,kq,lg). (63) 

We have the relations 

^iunorderedi^Jg^^g) = \(p3{iqJg^h) ^^3{iq^hJg)) ^ 

^iunorderedih^gA) = \{p3{h J g A) + HUg^h^ig) + ^iih^h^ J g)) ^ (64) 
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as well as 

lO 



^4,morderedi^q^ jg^^g^h) = 4 g^^g^h) ^^^ih^h^^g^ J g) ^^^ih^ J g^h^^g) 



^4,unorderedih^jq^^q^h) ~ ^4,qqqg{^qT Jqj^q^h) ^4,qgqq{^qjh^ Jq^^q)' 

The integrated antenna functions are denoted by a calligraphic letter, the generic ones are 

X^{sijk), xj{sijk), X^{sijki). (66) 
The integrated antenna functions are related to the unintegrated antenna functions by 

xi = S-'l/' £ fd<\>unresXi, ^4° = 5" £ f d<\>unres (67) 

orderings orderings 

Each antenna is given as a product of a pref actor P or (P, a symmetry factor S, a colour factor C 
and a kinematical factor K ox "K: 

X = PSCK, X = 'PSC^K. ■ (68) 
The prefactors are given for the unintegrated antenna functions by 

87rav, forX^ 

87ca,(|), forXj (69) 
(87ra,)% forXO 

For the integrated antenna functions the prefactor is given by 

r f , forxO 
'P = forx,' (70) 

The symmetry and colour factors are equal for the unintegrated and integrated antenna functions 
and are given in table 2. The normalisation of the kinematical factors K and ^ agrees with the 
antenna functions listed in ref . [37] . The introduction of the symmetry factor ensures for example 
that the antenna functions A3, D3 and or A4 and D'^ have the same soft singularities. As an 
example I list here the ^g^-antenna function in the integrated and unintegrated form: 

A%„jg,k,) = 8m,— f2^ + (l-8)^ + (l-8)^-28Y 



-e 



1 3 19 7 2 
? + 2^ + T-12" 
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+ 0(8). (71) 





AO DO eO fO go 


^3 ^3 ^3 ^3^ ^3 


aO 5O cO Z)0 


s 
c 


1 1 Ik 1 Ik 

Ca Ca 


1 i 1 i 1 

]^ ]^ Ik J Ik 

Ca Ca 


1 1 i i 1 
11111 



Table 2: The symmetry and colour factors for the antenna functions. 



The one-loop three-parton antenna functions (j, 7, k) are ultraviolet divergent and need to be 
renormalised. I will assume that these functions are all renormalised at a scale fi. The relation 
between the renormalised and bare one-loop three-parton antenna functions is given by 



and an identical relation holds between and ^K^hare- ^^^^ useful to know the infrared 
poles of the one-loop antenna functions. They are given by 

^3(^,7,^)|/M.. = {^3 - 2^3 (%■) -NfE^Sij) - 'D^{Sjk)-Nf'E^{Sjk) 

^^c ) IR-Poles 

D3{iJ,k)\ii,_p^i^^ - {2©0(.,;,) +2A^^2;0(5,,,) - ©0(5,,.) -Nfj:^isij) - 'D^Si^) -Nf^^{sik) 

-^,'>{sjk) - 2Nfg^{sjk) } dO(/, j, k) , 

+^^lisjk)]EliiJ,k) 



c 



(73) 

IR-Poles 



3.4 Phase space mappings 



The antenna functions in the subtraction terms are accompanied by amplitudes and observable 
functions, which are evaluated with configurations where one or more partons have been re- 
moved. The momenta of these reduced configurations are related by a map of momenta to the 
original configuration. The momentum map has to satisfy momentum conservation and the on- 
shell conditions. We first consider the case where a single particle is removed. Let [puPj.Pk) 
be a set of momenta corresponding to an antenna, such that particle j is emitted by the antenna 
formed by particles / and k and denote the mapped momenta of the reduced configuration by 
{PiiPk}- There are several possibilities to relate the mapped momenta {pupk) to the original 
momenta {pi,Pj,Pk)- One such possibility is [99]: 

. _ {^ + p)sijk-'2-rSjk ^ {l-p)sijk-2rsij 

^' 2{sijk-Sjk) ^' 2(sijk-Sij) 

« = (' -P>^^>-^(' -'■^^V . + d (74) 

2{sijk - s jk) 2[Sijk - Sij) 
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where 

p= ll+4r{l-r)^, r=— (75) 

V ^ijkSik Sij + Sjk 

This choice interpolates smoothly between the singular limits Pi\\pj and Pj\\pk- Setting however 
r = 1 yields the Catani-Seymour choice [48] 

Pi = Pi + Pj — Pk, Pk = — Pk, (76) 

^ijk ^ij ^ijk ^ij 

which can be used in the Pi\\pj limit, but not in the Pj\\pk limit. This mapping is symmetric in 
(i, j). On the other hand, setting r = yields 

Pi = — Pi , Pk = — Pi + Pj + Pk, (77) 

^ijk ^ jk ^ijk ^jk 

the Catani-Seymour mapping for the opposite case, which can be used in the pj \ \pk limit, but not 
in the Pi\\pj limit. This mapping is symmetric in {j,k). For the three-parton antenna functions, 
which have singularities in Pi\\pj and Pj\\pk we use the mapping eq. (74). This concerns the 
antenna functions 

kg), Fi{ig (78) 

The antenna functions and G3 have on the other hand only singularities for Pj\\pk and are 
symmetric in (7,^). We therefore use the mapping eq. (77), which respects the symmetry in 
(7,^) for the antenna functions 

E\{iqjq,kq), d^{igjq,kq). (79) 

For the double unresolved subtraction terms we have to relate four momenta {pi, Pj, Pk, Pi) to 
two hard momenta (pi.pi). The 4^2 mapping can be realised through two 3^2 mappings. 
There are two possibilities for four-parton antennas with an ordered singularity structure. The 
first possibility is to combine first {pi,Pj,Pk) into {pi,Pk) and then to combine {piiPkiPl) into 
{Pi,Pl)- 

{Pi,Pj,Pk,Pi) {Pi,Pk,Pi) {Pi,Pi)- (80) 

The second possibility is to combine first {pjiPkiPl) into {pj,Pi) and then to combine {pi,Pj,Pi) 
into {pi,pi): 

iPi,Pj,Pk,Pi) iPi,Pj,Pi) (Pi, Pi)- (81) 

As a criteria which mapping should be used we take that the smallest antenna should be recom- 
bined first. As a measure for the smallness of an antenna X^{iJ, k) we take the product 

SijSjk- (82) 
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In the limit where these invariants are small this quantity is proportional to the transverse mo- 
mentum of the emitted particle j. For the order {i,j,k,l) we therefore compare SijSji^ with sjj^si^i, 
or equivalently Sij with Ski- If Sij is smaller, one uses eq. (80), if on the other hand s^i is smaller 
one uses eq. (81). This procedure is used in combination with eq. (74) for the antenna functions 

^4{^qjjgjkg,lq), Al^^^(iq,jg,kg,lq), D^(iq, jg,kg,lg) . (83) 

For the antenna functions 

^Ai}qi jqi^qi^^i CAi}qi jqi^qi^q)i EA,qqqgi}qi Jqi^qi^g)- (^'^) 

we use eq. (76) or eq. (77) for the first mapping. For the antenna function 

E4,sci^qJq^kg^h) (^5) 

we always map {j,k,l) first with eq. (74) followed by a mapping {puPjiPi) with eq. (77). For 
the antenna function 

E4,qgqq(^T jg^^qJq) (^6) 

we map in the case Sij < s^i first («, j, k) with eq. (74) followed by a map (/?,■, Pk,Pi) with eq. (77). 
In the case Sjj > s^i the momenta {j^k^l) are mapped first with eq. (77) followed by a map 
{PuPjiPi) with eq. (74). These mappings basically follow the singularity structure of the antenna 
functions. 

A final remark on the notation for the arguments of the amplitudes: The particles in an 
amplitude are usually labelled by consecutive numbers. We have denoted the mapped momenta 
with a tilde, keeping the original numbering untouched. The mapped event has one or two 
partons less and leaves therefore a gap in the consecutive numbering. It is convenient to restore 
the consecutive numbering by relabelling all particles in the event. The momenta after this 
relabelling will be denoted with a prime. As an example consider the five-parton amplitude 

^f)igggq{PQ,PuP2,P3,PA) together with the mapping (pupi.Ps) ^ {puPs)- This yields 



which will be denoted by 



^4!lggq{P0,Pl,P3,P4) (87) 



41gqip'0^P'l^P2,P3)- (88) 



Therefore 



Po=Po, p'l^Ph P2^P3, p'a^P3- (89) 



The subtraction terms to the amplitude A ^l-, , (pQ ,pi,P2iP3,P4) proportional to correspond 

to a. g ^ qq splitting and we obtain as the reduced amplitude ^A^laan- convention we place 
the merged gluon immediately after the anti-quarks, if the spectator was an anti-quark. We place 
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the merged gluon immediately before the quarks, if the spectator was a quark. Therefore for 
£5(0,1,3) we have 

^4!qggqiP0,Pl,P2,P3) = ^i^lggqiPo, P3, P2, P4) , (90) 

while for £0(4, 3,1) we have 

^4!lggqiP0,P'uP2,P3) = ^4!qggq{P0,P2,PhP4)- (91) 

3.5 The hybrid method of subtraction and slicing 

Let us recall eq. (54) for the subtraction method at NNLO: The NNLO contribution is written as 
with 

(o(3))r- = /(of).of).of)o.ar.of)o.a^^)), 

^oi^)^NNLo ^ I ^^(3) ^^(0) _ ^(3) ^^^(0,1) _ of ) o Jaf . (93) 

The strongest conditions, which we can put on the form of the subtraction terms, are as follows: 

(1) The explicit poles in the dimensional regularisation parameter e of the three-parton contri- 
bution {o^^^)2^^^ cancel for each point of the three-particle phase space. 

(2a) The explicit poles in the dimensional regularisation parameter e of the four-parton contri- 
bution (^o^^^)^^^^ cancel for each point of the four-particle phase space. 

(2b) The phase space singularities of the four-parton contribution l^o^^^)^^'^^ cancel point-by- 
point in the four-particle phase space. 

(3) The phase space singularities of the five-parton contribution (^o^^^)^^^^ cancel point-by- 
point in the five-particle phase space. 

Due to spin-correlations in the coUinear limit the spin-averaged antenna function do not lead to 
a point-by-point cancellation of the phase space singularities. Therefore we have to relax the 
requirements on the subtraction terms. In addition we relax also the requirement on the cancella- 
tion of the explicit poles in the four-particle phase space. In a hybrid scheme of subtraction and 
slicing the last three items above are replaced by 

(2a') The explicit poles in the dimensional regularisation parameter e of the four-parton contri- 
bution ( 0*^3^)4^^^'^ cancel after integration over unresolved phase spaces for each point of 
the resolved phase space. 
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(2b') The phase space singularities of the four-parton contribution l^o^^^)^^^^ cancel after az- 
imuthal averaging has been performed. 

(3') The phase space singularities of the five-parton contribution l^o^^^)^^^*^ cancel after az- 
imuthal averaging has been performed. 

The azimuthal average is not performed in the Monte Carlo integration. Instead a slicing param- 
eter T| is introduced to regulate the phase space singularities related to spin-dependent terms. It is 
important to note that in this scheme no numerically large contributions proportional to a power 
of Inri are cancelled between (o^^))^^^^^ ^q{3)^nnlo ^q{3)^nnlo_ Each contribution itself is 
independent of T] in the limit T] — 0. 

Let me discuss the cancellation of the explicit poles in the four-particle phase space (item 
2a') in more detail: We consider contributions of the form 



(94) 



The functions <y^{i,j,k, ...) contain the explicit poles in the dimensional regularisation parameter 
£. The dots indicate that these functions may depend on partons other than partons j or k. Item 
2a' states that we do not require that the explicit poles cancel point by point in the four-parton 
phase space, which is given by J(|)4 = d^s^-id^s^unres^ijk- Instead we only require that these poles 
vanish after integration over the unresolved phase-space d^unresjjk- This is sufficient to ensure 
that the expression in the square bracket is finite and the integration over the three-parton phase 

(3) 

space d^2 can be performed in four dimensions. The observable O3 ^ and the matrix element 

I Jlj^'* P can also be evaluated in four dimensions. 

Now we can distinguish two important cases. In the first case the expression in the square 
bracket vanishes after integration over the Z)-dimensional unresolved phase space d^unres,ijk to all 
order in e. This case is relevant for the subtraction terms in the main part of this paper. Obviously, 
since the square bracket vanishes to all orders in e we have in particular that the term of order £^ 
vanishes. Therefore we do not obtain any finite terms in this case. 

In the second case only the pole terms vanish after integration over the D-dimensional unre- 
solved phase space d^unres,ijk- In this case eq. (94) will lead to finite terms. This case is discussed 
in appendix A and only relevant to the alternative subtraction terms presented there. 

4 The subtraction terms 

4.1 The tree-level single unresolved subtraction terms 

We first give the NLO subtraction terms Jttg"*'''^ and da^"^^^ . The first one gives the complete 
subtraction term for the NLO coefficient, the second one enters the NNLO coefficient. We have 



^<S';. - [dO(0, 1,2) +D«(0,2,1) +2)0(3, 1,2) +2)0(3,2,1)] (95) 
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INr 



1-[aO(0,1,3)+aO(0,2,3)]1o 



SI 



(0) 



J(|)4, 



3,qq'q'q 



7 + ^1 V [^3 (0, 1,2) +£?(1,0, 3) +£?(2, 3,0) +E?(3,2, 1)] 



IN, 



+ -^[E^{0,l,3)+El{l,0,2) + E^{3,2,0)+El{2,3,l)] 



SI 



(0) 



(i(j)4. 



The amplitude Sif^ is evaluated with momenta obtained from eq. (74), which depend on the 
partons, which form the antenna. Integration over an one-parton phase space yields 



/ 



dec 



single 



'^^3,qggq^^^3,qq' 



>le \ 
^q'q ) 



(96) 



N. 



' [2)30(,0i) +A^/£0(5oi) + 2)3^^12) +%£3 (^12)] -2^-^3(^02) 



SI 



(0) 



#3. 



For the subtraction term doC'^'^^^'^ we have to take into account colour correlations. For the sub- 
process e+e qgggq we first define a two-component vector in colour space obtained from the 
colour-ordered amplitudes of e^e" — > qggq: 



^4\q0,gl,g2,q3) = eVco 



2„2„_ ( Af\qo,gi,g2,q3) 



^4°n?0,g2,gl,?3) 



(97) 



For the subtraction term we have 

1 single _ 1 A^c ~ 1 
""4,9ggg9 3 8 

^0 



[Z)0(0, 1,2) +Z)0(0,2, 1) +^0(0,3, 1)] o.^^' 



mfN^-i -1 



+ [dO(0, 2, 3) +D^(0, 3,2) +dO(0, 1, 3)] o A, 



(o)t 



-1 -1 

1 -1 
-1 N}-1 



SI 



(0) 



SA 



(0) 



+ [F0(l,2,3)+F«(l,3,2)+F«(2,l,3)]o5ft|^^' ) 
+ [D§(4,2, 1) +D§(4, 3, 1) +i)?(4, 1,2)] o 5f ^ 



-1 -1 
-1 N}~\ 



5f 



+ [D0(4,3,2)+D§(4,l,3)+D?(4,2,3)]o5i' 



-;(0)t / A^c -1 -1 \ ;;(0) 



-1 -1 



S^A 



+ 



[A§(0,l,4)+A§(0,2,4)+A«(0,3,4)]o5fM ^ 1+- 



S#5. (98) 



Af2 Af2 



For the sub-process e^e qqq'q'g we have two subtractions terms, one involving the two- 
quark two-gluon amplitudes of e~^e~ — > qggq, the other one involving the four-quark amplitudes 
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of e'^e — > qqq'q' 



doc 



single i/oc*'"^'^ 



,+da: 



single 



'4,qq'gq'q-^qggq ^ 4,qq'gq'q^qq'q'q- 



(99) 



We have 



4,qq'gq'q^qggq 



1 Nf-l 



1 



(£»(0, 1,3) +£0(1,0,4) +£30(3,4,0)+£0(4,3,l)) 



+2^ (E?(0, 1,4) + £0(1,0,3) +£?(4,3,0) + e;(3,4,1)) 

+ 2(1 + ^ 
'4 2 



-1 A^,^-l 



5f 



(G§(2, 1 , 3) + G§(2, 0,4)) + 2I (G§(2, 0, 3) + G§(2, 1 , 4)) 



(100) 



The asymmetry in the colour factors for the i? -terms is just related to way the partons are labelled 
in the amplitudes. For the subtraction term da''!"^i^,^,^--,,, we introduce a two-component 
vector in colour-space obtained from the colour-ordered amplitudes of e+e qqq'q': 



-5(0)/- _/ / N 2 2 



%r(^0,^l,?2'?3) 

^fiavxTiq'iAQA'i^q^) 



(101) 



We then have 



dec 



single 



'A.qq'gq'q^qq'q'q 



[A§(0,2,3)+A§(l,2,4)]oxJ 



(o)t 



1 



+ [aO(0,2,4)+aO(1,2,3)]oxJ 

+ [A0(0,2,l)+A0(3,2,4)]oxr/ ^ 
Integration of these subtraction terms yields 



1 

1 



(0)t ( Nc 

_ 1 _ J. 



7(0) 



(102) 



1 + 



1 



1 +^ ^ 



, single ,j^single \ _ 

^^4,qgggq '^^4,qq'gq'q^qggq J " 



[©30(501) + 2)30(523) + A^/£3°(*0l) +Nf'E^{s23)] K 
+ [©30(502) + 2)30(513) + A^/2;3'(*02) +%2;3'(*13)] 
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(o)t / A^c -1 -M :;(o) 



-1 -1 
-1 -1 

-1 N^-l 



(0) 



1 



+ [^0(^03)] ^ I _ >^ ) if ^ #4, (103) 



2 1 



+ [;.?(^0.) + ^3°(^23)] 2f * ( 1 f ^ ' ) \ d^4. 



In the first equation terms proportional to £3 have been equally distributed between the two 
colour matrices. 

4.2 The one-loop single unresolved subtraction terms 

For the subtraction term do!^"^ we have two contributions corresponding to the sub-processes 
e+e^ qggq and e+e^ — > qqq'q' . Each of these contributions will have terms proportional 
to the one-loop interference term 2Re ^f'^'^ R^^"^ and the tree-level antennas as well as terms 

proportional to the Born matrix element |-^^3'^'' P and one-loop antenna functions X3 . We therefore 
write 

= dai^'P+dal'Z^^u- (104) 

3,qq'q'q 3,qq'q'q,a 3,qq'q'q,b ^ ' 

We have 

= ^{y W(0,1,2)+d0(0,2,1)+d0(3,1,2)+d0(3,2,1)] (105) 

[A§(0,l,3)+A§(0,2,3)]}o (^f t^3W + ^3Wt^^^^ 
= ^{^[^3(0,1,2)+Di(0,2,l)+Di(3,l,2)+Di(3,2,l)] 



-i[Ai(0,l,3) + (A^2^l)Ai,,,(0,l,3) 

2 



+Ai (0, 2, 3) + (Ar2 + 1) A^,,, (0, 2, 3)] } o 



#4, 
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^<Iy,,« = {Q + ^)|[£«(0,l,2)+£?(l,0,3)+£3«(2,3,0)+£3«(3,2,l)] 
+1 V [^3 (0, 1,3) +£?(1,0,2) +^0(3,2,0) +E?(2,3, 1)] 



d^tm',,b = { + ^) T (0, 1 , 2) + £3' ( 1 , 0, 3) + E\ (2, 3,0)+ (3,2,1 )] 

+ ^^[£i(0,l,3)+£i(l,0,2)+£3i(3,2,0)+£3i(2,3,l)]|o|;lf |'j^4. 
Integration over an one-parton phase space yields 

1 1 

^ [2>3^^0l)+A^/%^^0l) + 2>3^^12)+%£3M] -2^-^3(^02) 
1^ [2)3l(,Ol)+A^^£3l(,o0 + ©3n^l2)+A^/£3H^12)] [^3'W + (A^c' + ^ 

#3- (106) 



+ 



•^3 



(0) 2 



4.3 The iterated single unresolved subtraction terms 

Whereas the subtraction term da!^"^ approximates the one-loop matrix element in singular phase 
space regions, there are on the four-parton phase space additional terms, descending from the 
five-parton phase space: 

Jaf ^'^ + da'^^"""' + dall"''''''^ + Jaf (107) 

These terms themselves are singular in single unresolved phase space regions, except for the 
term da''^'^\ which is integrable in all limits. The term da^^'^'^'*^ is a subtraction term for these 
contributions. In the unintegrated version da^"'^"^* lives on the four-parton phase space and 
involves one unintegrated three-parton tree-level antenna function together with one integrated 
three-parton tree-level antenna function, which however does only depend on the hard three- 
parton kinematics. Therefore this integrated antenna function factorises from the integration 
over the unresolved phase space and the integrated version is just a product of two integrated 
three-parton tree-level antenna functions with three-parton kinematics. We have 

= ^ {t 1'2) +Z)0(0,2, 1) -Z)0(3, 1,2) -^0(3,2, 1)] (108) 
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m 1 1 



4 2 4 



+ - [dO(0, 1,2) +Z)0(0,2, 1) +D«(3, 1,2) +^^(3,2, 1)] oji^{so,2') 



+^ [aO(0, 1,3) +A0(0,2, 3)] o ^2)0(50.10 + ^£3°(^0'1') + \^!{si'2') + f £3°(^1'2')) 
+ + i) ^'^^ +^3(0,2,3)] o^K.o^^,)} o 



^ + \^)[^3(0,1,2)+^?(1,0,3)-e0(2,3,0)-£3«(3,2,1)] 



+i [^0(0, 1,3) +^0(1,0,2) -£0(3,2,0) -^0(2,3, 1)] j o (^—^'.{so^v) + ^2?3°(^i'2')) 

[E^{0,l,2)+E^{l,0,^)+E^{2,3,0)+E^{3,2,l)] 



— - + - 
4 2 4 



1 I 

4 2 



+i [£0(0, 1,3) +£«(1,0,2) +£0(3,2,0) +i?3«(2,3,l)] 



Integration over the one-parton phase space yields 



O-*?3(*0'2') \° 



(0) 



J{|)4. 



/ 



dec 



product 



-\^'{s,,f - U',{snf + 2?30(50l)2>3°(^12) + U'^is^if 



N^Nf 



[-2>?(5oi)£?(^oi) - 2>?(^i2)£?(^i2) + ^!isoi)'El{sn) + 2)?(^i2)2:?(*oi)] 



+ ^ [-'D3M^3M - 2)3°(512)^3 (^02) +^3'(*02)'] 
+ ^ [-£3'(^0l)^3°(^02)-£30(.12)^3'M j 



J(^3. 



(109) 



Remark: The term 



[22;30(.oi)2;3'(^12) - £3°(^0l)' - £3' (^12)'] 



J? 



(0) 



ii(l)3 



(110) 



is finite. Therefore eq. (109) can be re- written in a form, such that the antenna function £3 occurs 
only in the combination D^+Nf'E^ modulo finite terms of the form (110). 

4.4 The colour-connected double unresolved subtraction terms 

We consider now the subtraction terms for the double real emission contribution. We start 
with the subtraction terms involving the four-parton antenna functions. For the subtraction 
term d(X'^°^^^^ we have two contributions corresponding to the sub-processes e^e" — > qgggq 
and e~^e~ — > qqq'q'g. We have 
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da: 



double 



1 

6 



£ ^[Dl{Q,iJ,k)+Dl{A,iJ,k)\ 



J,j,ke{l,2,3},i^j^k 



+ I X [-4.c(0,^7,4)] + i [-A0(0,/,7,4) -A0,,(0,/,7,4)] 
Ue{l,2,3},iW ^ 



-^Kc(0,/,7,4)] 



(0) 



#5, 



(111) 



j^double 



1 Nf-l 
- + — 

4 2 



(112) 



4 



1 1 



4iV.2 



<,,,,(0,l,3,2)+4^^^^(l,0,4,2)+4^^^^(4,3,l,2)+£4%^(3,4,0,2) 

<,g,,(0,2,3, 1) +4,^^,^^(1,2,4,0) +4^^^^(4,2, 1,3) +£4^^^(3,2,0, 4)) 
i(4,,(0,l,2,3)+4jl,0,2,4)+4j4,3,2,l)+4,,(3,4,2,0)) 



^(^0(0, 1,3,4) +50(1,0,4,3)) 



+ 

+ 



1 
4 

4 



4 V<,g,g(0: 1,4,2) +^^^^^^^(1,0,3,2) +£4^^^(4,3,0,2) +£4"^^^^(3,4, 1,2) 
4,,,,(0,2,4, 1) +£4^^^(1,2,3,0) +£4^^^(4,2,0,3) +<,,,,(3,2, 1,4)) 



+ - 



1 1 



4iV.2 



(£4«,,(0,1,2,4)+£0,,(1,0,2,3)+£4V(4,3,2,0)+£4«,,(3,4,2,1)) 



4K 



^(fi«(0, 1,4,3) +50(1,0,3,4)) 



-\ (1 - [C5(0, 1,3,4) +C5(0, 1,4,3) +(^(1,0,4,3) + C5(1, 0,3,4)] | 



(0) 



#5- 



Integration over a two-parton phase space yields 



da' 



double 



7 double _i_jfy> 



^double 



^ W(^Ol) + 2)0(^12) - <,(.02)] +^ [£4' (^01) + £4°(M] 
+ \ [-^4(^02) -<,c(*02) -4^0(^02)] + ^ -!S4(*02) - ^^^tcM " ^£4^(^12) 



W„c(^02)+4ci'(502)] 



(0) 



#3- 



(113) 
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4.5 The almost colour-connected double unresolved subtraction terms 



For the subtraction term daf"^^' we have two contributions corresponding to the sub-processes 



e ' e 



qgggq and e'^e — > qqq q g. For the sub-process e'^e — > qgggq we have 



doc 



almost 



1 fAgl 

4 2 



[(-D^(0,1,3)-Z)0(0,2,3)-D0(0,3,2)+Z)0(4,3,2)+D^(4,1,3)+Z)0(4,2,3)) 
o(D«(0M',20-Z)0(3M',2')) 

+ (-D0(4,2,1)-D0(4,3,1)-D0(4,1,2)+Z)0(0,1,2)+d0(0,2,1)+Z)0(0,3,1)) 

o(Z)^(3',2MO-dO(0',2MO) 

-(A«(0,1,4)+A«(0,2,4)+a0(0,3,4)) 

o (dO(0', 1',2') +Z)0(3',2', 10 +Z)^(0',2', l') +Z)«(3', l',20) 

+ (dO(0, 1,3) +D^(0,2,3) +Z)^(0,3,2) +Z)0(4,3,2) +^0(4, 1,3) +^0(4,2,3)) 

oAO(OM',30 

+ (d0(0,1,2)+d0(0,2,1)+d0(0,3,1)+d0(4,2,1)+d0(4,3,1)+d0(4,1,2)) 

oA«(0',2',3') 

+ (aO(0, 1,4) +a0(0,2,4) +a0(0,3,4)) o (aO(0', 1',3') +a0(0',2',3'))] 
+i[(A§(0,l,4)+A§(0,2,4)+A§(0,3,4))o(A§(0',l',3')+A§(0',2',3') 



-D0(0', 1',2') -Z)0(3',2', 1') -Z)«(0',2', l') -£>^(3', l',2'))] } o 



(0) 



#5- 



(114) 



For the sub-process e+e qqq'q'g we have two subtractions terms, corresponding to an inter- 
mediate qggq or qqq'q' state: 



jf^almost 
^^3,qq'gq'q 



jf^almost I jfyCilmost 

^M'gq'q^qggq^m ^M'g<iq^m'q'q^m- 



We have 



almost 

3,qq'gq'q^qggq^qgq 



lAgl 
2 4 2 



(115) 



(116) 



^ + \^) (-^3>,l,3)-^3°(l,0,4)+£30(4,l,3) + £30(3,0,4)) 

+^(-£30(0,l,4)-£30(l,0,3)+£0(3,l,4) + £0(4,0,3)) 

o(D«(0',l',20-Z)0(3M',20) 
1 %-l 



+ 



4 + - 2 



(-£30(4,3, 1) -£3^(3,4,0) +£30(0,3, 1) +£30(1,4,0)) 



+1 (-£0(3,4,1) -£30(4,3,0)+£30(0,4,l)+£0(l,3,0)) 
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o(D0(3',2',l')-i)0(0',2',l')) 



+ 



1 Nf-l 
4 2 



{El{0,l,3)+El{l,0A)+E^{4,l,3)+E^{3,0A)) 



+i (^0(0, 1,4) +i?3"(l,0,3) +i?«(3, 1,4) +i?3°(4,0,3)) 



oA0(0',1',3') 



+ 



1 I Nf-l 
4 2 



(£§(0,3,1) + £;(1,4,0)+E?(4,3,1) + E?(3,4,0)) 

oA0(0',2',30 



+^{E^{0A,l)+E^{l,3,0) + E^{3,4,l)+E^{4,3,0)) 



(0) 



#5- 



da: 



almost 

3,qq'gq'q-^qq'q'q-*qgq 



1 A^f-1 
- + — 

4 2 



A«(l,2,4)-lA0(0,2,3)-iA0(3,2,4)-iA0(0,2,4))o£30(0',l',2') 



+ (^A0(l,2,4)-lA0(0,2,3)-iA0(0,2,l)-iA«(l,2,3))o£0^^^^ 

+ (^A0(0,2,3)-lA0(l,2,4)-iA0(3,2,4)-iA«(l,2,3))o£0(l',0',3') 



+ (^a0(0,2,3) - ^a0(1,2,4) - iA0(0,2, 1) - ^a0(0,2,4)^ o£0(3',2', l') 
+i(^(^A0(l,2,3)-iA0(0,2,4)-iA0(3,2,4)-iA0(0,2,3))o£3>^1^3^^ 
+ (a0(1,2,3) - ^a0(0,2,4) - 1a0(0,2, 1) - ^a0(1,2,4)) o£30(3',2',0') 



+ (A0(0,2,4)-lA0(l,2,3)-iA0(3,2,4)-iA0(l,2,4))o£30(l',0',20 



+ ( A«(0,2,4) - ^a0(1,2,3) - iA^(0,2, 1) - ^A«(0,2,3) ) o£30(2',3', l'^ 



1 

+4 



1 Nf-l 

- + — 

4 2 



(2(A§(0,2,1)-A§(0,2,3)+A§(3,2,4)-A§(1,2,4)) 

o(£30(0',l',20+i?3^l',0',30+£3°(2',3',0')+£3°(3',2',l')) 
-A«(0,2,4) o (£0(0', l',20 +£3«(3',2', l')) 
-A0(1,2,3)o(£30(1',0',3')+£3°(2',3',0'))) 

+^(2(a0(0,2,1)-A0(0,2,4)+a0(3,2,4)-A0(1,2,3)) 

o(£30(0',1',3')+£3°(1',0',2')+£3«(3',2',0')+£3°(2',3',1')) 
-A0(0,2,3) o (£30(0', l',3') +£3«(2',3', l')) 
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-A0(1,2,4)o(£0(1',0',20+£3°(3',2',0')))]} 



(0) 



These subtraction terms are integrated over an one-parton phase space. We obtain 



almost 

3,qq'gq'q-^qggq-*qgq 



(117) 



Ni 1 



4 2 



[(D§(0, 1,2) -D0(3, 1,2)) (2>?(523) + A^/2;?(^23) - ©?(^02) -%2;?(^02)) 

+ (dO(3,2, 1) -D^(0,2, 1)) (2)«(50l) +A^/£3°(50l) - ©3«(5i3) -%£3"(*13)) 

- (dO(0, 1, 2) +d0(3,2, 1) +dO(0, 2, 1) +^0(3, 1,2)) ^^(^os) 

+A0(0, 1,3) (©0(502) +iV/£3'(^02) +2^3'fe) +%£3°(^23)) 

+a0(0,2, 3) (©0(501) +A^/%«(5oi) + ©3'(^i3) + A^/£3 (^13)) 
+ (a0(0,1,3)+a0(0,2,3));130(503)] 

+1[(a0(0,1,3)+a0(0,2,3)-D§(0,1,2)-D0(3,2,1)-D0(0,2,1)-D§(3,1,2)) 

2 



■^3(^03)]}o 



(0) 



/ 



j^almost 



N 



1 Nf-l 

- + — 

4 2 



£0(0, 1,2) f;l30(.i3) - ^^0(^02) - ^^3^(^23) 



+£0(1,0,3) f^30(.o2) - U'sM - UU^23) - U'si^n) 



+£30(3,2, 1) ( ^30(502) - U!{si3) - Uliso^) - UU^03) 



+ 1 (eI{0, 1,3) (j^^{si2) - i^3°(^03) - U^{S23) - U!{S02) 



+£30(3,2,0) (jl',{sn) - U'3M - U^M - U'3M 



+£30(1,0,2) f^30(.03) - ^-^3°(^12) - ^-^3'(^23) " ^-^3°(^13) 



+£30(2, 3, 1) [^30(503) - ^^3"(^12) - ^-^3" (^01) - ^^3°(^02) 



+ 



^+^^](2(£30(0,l,2)+£30(l,0,3) + £30(2,3,0)+£30(3,2,l)) 



4 2 

(^30(.0l) -^3°(^02) +^3°(^23) -^3'(^13)) 

-(£30(0,l,2)+£30(3,2,l))jl30(503)-(£30(l,0,3)+£30(2,3,0))^30(.i2)) 
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+^ (2 (£0(0, 1,3) +£3"(1, 0,2) +£0(3,2,0) +£0(2,3, 1)) 

(^30(.0l) -^3°(^03) +^3°(^23) -^3'(^12)) 

-(£?(0,1,3) + £30(2,3,1))^?(502)-(£?(1,0,2) + £;(3,2,0))^30(5i3))]}o 



(0) 



4.6 The iterated double unresolved subtraction terms 

For the subtraction term da'^^'''^^'^ we have two contributions corresponding to the sub-processes 
e'^e~ — > qgggq and e'^e~ — > qqq'q'g. For the sub-process e'^e~ — > qgggq we have 



da 



iterated 

3, qgggq 



N: 



(118) 



^[(d0(0,1,2)+d0(0,2,1)+d0(0,3,1)+d0(4,3,2)+d0(4,1,3)+d0(4,2,3) 

+F30(1,2,3)+F30(1,3,2)+£30(2,1,3))o(d0(0',1',2')+Z)^(3',2',1')) 

+ (d0(0,2,3)+d0(0,3,2)+d0(0,1,3)+d0(4,2,1)+d0(4,3,1)+d0(4,1,2) 

+F30(1,2,3)+F30(1,3,2) + F30(2,1,3))o(d0(0',2',1')+Z)^(3',1',2'))] 

+\ [(^ + (^'(0, 1,4) +a0(0,2,4) +a0(0,3,4)) o (a0(0', 1',30 +a0(0',2',3')) 

-(a0(0,1,4)+a0(0,2,4)+a0(0,3,4)) 

o (d0(0', 1', 2') +Z)^(3',2', 1') +Z)0(0',2', l') +Z)0(3', l', 2')) 

-(dO(0,1,2)+D?(0,2,1)+dO(0,3,1)+dO(0,2,3)+dO(0,3,2)+D?(0,1,3) 

+d0(4,2,1)+d0(4,3,1)+d0(4,1,2)+d0(4,3,2)+d0(4,1,3)+d0(4,2,3)) 

o (a0(0M',3') +aO(O',2',30)] } o |^3(«)|' J(^5. 

For the sub-process e^e~ qqq'q'g we have two subtractions terms, corresponding to an inter- 
mediate qggq or qqq'q' state: 



jf^iterated 
^^3,qq'gq'q 



J ^iterated _ _ _|_ ^/v" 

3,qq'gq'q-^qggq^qgq '^^3,qq'gq'q-*qq'q'q-*qgq- 



^iterated 



(119) 



We have 



Jterated 

3,qq'gq'q-*qggq-*qgq 



1 - 1 



(£30(0, 1,3) +£?(1,0,4) +£?(3,4,0) + £0(4, 3, 1)) 



+i (£30(0, 1,4) +£30(1,0,3) +£30(4,3,0) +£30(3,4, 1)) 

(d0(0',2', 1') +dO(3', 1',2')) - I (aO(0', 1',3') +aO(0',2',3')) 



+- 



4 



1 Nf-l 

- + — 

4 2 



(g0(2, 1,3) +g0(2,0,4)) + 1 (g0(2,0,3) + g0(2, 1,4)) 
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o (Z)0(0', l',2') +D^(3',2', 1') +d0(0',2', 1') +£>^(3', l',2'))}o 



(0) 



(120) 



jf^iterated 



1 

- + ^^ 

4 2 



(A»(0,2,3)+A^(l,2,4)) 



+^(-A?(0,2,4)-A?(1,2,3)-2A0(0,2,3)-2A§(1,2,4) + 2A§(0,2,1) + 2A?(3,2,4)^ 
o [£0(0', l',20 +£30(1',0',3') +£3°(2',3',O0 +£3^3',2', l')] 



+- 



(A§(0,2,4)+a0(1,2,3)) 



+-(-A0(0,2,3)-A0(1,2,4)-2a0(0,2,4)-2a0(1,2,3) + 2a0(0,2,1) + 2a0(3,2,4)) 



o [£0(0', l',3') +£30(1',0',2') +£0(3',2',0') +£3^2^,3', l')] } o 



SI 



(0) 



#5- 



(121) 



These subtraction terms are integrated over an one-parton phase space. We obtain 



jf^iterated , jf^iterated 



(122) 



{'D^isoi) +Nf'E^{soi) + ©30(523) +Nf'E^{s23) + !F3\sn)+2Nfgl{sn)) 
+ (d0(0,2,1)+d0(3,1,2)) 

(©30(502) +%£3°(*02) + ^!isi3) +Nf'E^isi3) + !Fi\sn)+2Nfglisi2))] 

+i [- (dO(0, 1,2) +dO(3,2, 1) +dO(0,2, 1) +dO(3, 1,2)) ^30(503) 

- (aO(0, 1,3) +aO(0, 2,3)) (©3^(501) +Nf'E^{soi) + 'D^{so2) +Nf'E^{so2) 

+ ©30(513) +iV^2;30(5i3) + ©3°(*23) + A^/'E3°(*23)) 



+ 1 + 



i^) W(0,1,3)+a0(0,2,3))^30(503) 



(0) 



I 



dec 



iterated 

■3,qq'gq'q-^qq'q'q-^qgq 



1 A^f-1 
- + — 

4 2 



[£30(0, 1,2) + £30(1,0,3) +£30(2,3,0) + £30(3,2,1)] 



2 

^(i?30(502)+^3°(^13)) 



+ i (-.^30(503) -^3°(*12) -2^?(*02) -2.;??(5i3) +2^30(501) +2.^30(523)) 

+i [£30(0, 1,3) +£30(1,0,2) +£30(3,2,0) +£30(2,3, 1)] 

^(;?30(503)+^3°(^12)) 
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(0) 



d(^4. 



4.7 The soft subtraction terms 

In the four-parton channel the poles from ^04^^ and da^^"^''^ cancel each other. The other terms 
which we defined up to now in the four-parton channel are 

da^^'P+daP""^""-! daf + J d^ 



iterated 



(123) 



We now study the infrared pole structure of the combination of these terms. In the four-parton 
channel there are two sub-channels, the qggq sub-channel and the qq'q'q sub-channel. We first 
consider the qggq sub-channel. For the infrared poles proportional to the unintegrated three- 
parton antenna function 1)3(0, 1 , 2) we find 



Po/e5,D?(0,l,2) 



2 4 2 



Z)^(0,l,2) (124) 



[2^3(^02) - 2>3 (^02) - ^3°(%) + ^3^*03) - 2)3^^23) + 2?3 (*23)] « 



(0) 



d(^4 



Poles 



There are three similar terms with poles proportional to 1)3(0,2, 1), 1)3(3,2, 1) and /)3(3, 1,2). 
For the infrared poles proportional to the unintegrated three-parton antenna function A3(0, 1,3) 
we find 



Poles A{Q,\,'i) 



1 fNl 1 1 

2 VT2 



+ ^)a?(0,1,3) 



[--^3 (^03) + -^3 (*03) + ©3^(^02) - 2^3 (*02) + 2>3°(^23) ' 2^3(^23)} o 



(0) 



d^4 



Poles 

(125) 



There is a similar term proportional to A3(0, 2, 3). Let us now consider the sub-channel qq'q'q. 
For the infrared poles proportional to the unintegrated three-parton antenna function £3 (0, 1,2) 
we find 



4 2 



1 1 

[2^3 (^02) - 2-^?3°(*0l) - ^3^%) + -^3 (*02) + Sl^{sQ3) - ©3^*23) + ^3 (*23)] 

2 



Poles, E^(Q,\,2) 



' '^'-'^^?'e?(0,1,2) 



4 2 



SI 



(0) 



J(|)z| 



(126) 



Poles 
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There are three similar terms with poles proportional to £'3(1,0,3), M'(2,3,0) and £3(3,2, 1). 
In addition there are four terms with proportional to £3(0,1,3), £3(1,0,2), £3(3,2,0) and 
£3(2,3,1) for which one replaces 1/4+ {Nf — l)/2 by 1/4. We can re-write the terms in 
eq. (126) as 

2>?(*012) - 2Jll{soi) - J?3°(5o3) + -^3 (*02) + ^3 (^Ob) - 2^3 (^23) + -^3 (*23) = 
[2>?(*012) - - + ^?(*03) - 2>?(*23) + "^3 (*23)] 

+2[j?30(502)-^3°(^0l)]. (127) 



Terms with 



|£30(0,1,2) [^0(502) -^3'(^0l)]o 



(0) 



d^4 



(128) 



vanish after integration over the unresolved phase space for a momentum mapping symmetric 
in {pi,p2). On the other hand, the first term on the right hand side of eq. (127) has the same 
structure as the terms in eq. (124) and eq. (125). The integrated antenna functions Sl^is) and 
2)3 (5) have the form 



a, 

71 



1 C 

8^ e 



+ 0(8). 



(129) 



Due to the structure of alternating signs in eq.(124), (125) and (127) the coefficient c of the single 
pole in eq. (129) will drop out and only the double pole contributes. The double pole is related 
to soft gluons. Generically we have to consider terms of the form 



X^{a,i,b) \s^{sgi,)-S^{sab) -S^isaj) +sUsaj) -5?(^5,-) +S^{sbj) 



(0) 



#4 (130) 



with 



^3(^) 



OC.S 



71 \//^ 

Neglecting finite terms this is equivalent to 



■ + - + 6- 



12 



-71 



+ 0(£). 



(131) 



^X^{a,i,b)- 
71 '8 



,(0) 



In 



-In 



Sab 



SI 



(0) 



(132) 



The three-parton amplitude is evaluated with the momenta pa, pb and pj. The three-parton 
amplitude and the observable are independent of particle i. We now investigate the behaviour of 
eq. (132) if we integrate over the phase space of particle i. It is convenient to choose a specific 
frame, which we take as the centre-of-mass frame of pa + pb with pa and pa along the positive 
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z-axis and pj in the x — z plane. In this frame we have 



1 

2" 
1 



Pa = ;^V^(1, 0,0,1), 



Pb = -v/S^(l, 0,0,-1), 

Pa = £«(1, 0,0,1), 

Ph = £'/,(l,sin6/,cos(|),sin0fesin(l),cos6i,), 

Pj = Ej{l,smej,0,cosej). (133) 



One finds 



ln!^-ln^^ = \n( (1 + cosej)(l -cosOz,) y ^^^^^ 

5 Sab V 2 ( 1 — COS 6^, COS By — sin 6^ sin Qj cos ^) J ' 

We focus on the integral over the azimuthal angle ^ of particle / in this frame. In eq. (132) only 
the square bracket depends on the azimuthal angle (|). With the help of eq. (134) the relevant 
integral is given by 



271 



ln( .,. ''+''°''y''rr^''l itI, (135) 



2% J V 2(1 — cos6i,cos0, — sinG^sinG/cosd)) 




The integral equals 



J ^ \n( ^~ ^ ^^"^ ~ ^ (136) 

\ 1 — COS0f,COS0y + I COS0y — COS0t| / 

The integral is zero for Qj < 0^ but non-zero for Qj > 0^. Therefore the 1/e-poles of eq. (130) 
vanish after integration over the azimuthal angle for Qj < 0^. However, they do not vanish for 

Qj > 0^ and an additional subtraction term da^^^* is needed to cancel the explicit poles of 

da^^^'P + daP'"^"'' - J daf + J daf'''''^. (137) 

1 1 

We introduce the notation 

S%,j,k;a,b) = 87ca,^^, (138) 

together with the convention that the phase space mapping ( j, _/, k) — > (/, k) should be used. In the 
qgggq-chaaael the soft subtraction term can be taken as 



36 



4 2 



£ (S0(0, /, 0", l") - 8^(0, i, j; O", l") - S',(Q, i, j; l", l") 

L0j)e{(l,3),(2,3),(3,2)} 

-50(0, O', l') + ^^(O, /, j; O', 3') + ^^(O, i, j; 2', 3')) o ^^(O', l', 2') 

+ ('^3(0, ^ 7;0", 1") - 5«(0, ;;0^ 2'') - 5«(0, l", l") 

(U)e{(l,2),(2,l),(3,l)} 

-S^(0,/,7;0',1') + S^(0,/,7;0',3')+5^(0,/,7;1',3'))oD0(0',2',1') 

+ £ (50(/,7,4;l",2") -5§(/,j,4;0", l") -5?(/,7,4;0",2") 

(U)e{(l,2),(2,l),(l,3)} 

-5S(/,7^4;1^30 + 5§(/,7^4;0^l0+5§(/,7,4;0^3'))oD0(3^2^l') 

+ £ (s0(/,7,4;l",2")-50(/,7,4;0M")-50(/,7,4;0",2") 

Oj)e{(2,3),(3,l),(3,2)} 

-50(/,7^4;2^30+50(/,7^4;0^2')+50(/,7•,4;0^3'))oD0(3^1^2')] 

X^ + i) \ I (5^(0,/,4;0^2'0-5°(0,/,4;OM'0-S«(0,/,4;1^2") 

^ ^ U{1,2,3} 

-5§(0,/,4;0',3')+5§(0,/,4;0',2') + 5?(0,/,4;2',3'))oA?(0',l',3') 

+ £ (5?(0,/,4;0",2")-55(0,/,4;0",l")->^§(0,/,4;l",2") 

ie{ 1,2,3} 



-5^(0,/,4;0',3')+5^(0,/,4;0',1') + 5^(0,/,4;1',30)oA0(0',2',3')]} 
In the ^^'g^'^-channel we have 



(0) 



^mgqq 4 2 
1 

- + ^^ 

4 2 



(140) 



^ [(5^(0,2,3;0",n-5^(0,2,3;0",2")-5«(0,2,3;l",2") 



-5^(0, 2, 3; 0', 2') + 5^(0, 2, 3; O', 3') + 5^(0, 2, 3; 2', 3')) o £30(0', l', 2') 
+ (5^(1, 2, 4;0",1") -5^(1, 2,4;O^2'0 -5^(1,2,4; 1^2") 
-5§(1,2,4;1',3') + 5^(1,2,4; l',2') + S^(l,2,4;2',3'))o£30(l',0',3') 
+ (5^(0,2,3; 1",2") -5^(0, 2,3; 0", l") - 5^(0, 2, 3; O", 2") 
-5S(0, 2, 3; 0', 2') + 5§(0, 2, 3; O', l') + ^5(0, 2, 3; l', 2')) o E^{2' , 3', O') 
+ (5§(1,2,4; 1",2") -5§(1,2,4;0", 1") ->^§(1,2,4;0",2") 
-S0(l,2,4;l',3') + 50(1,2,4;0',1') + S0(l,2,4;0',3'))o£30(3',2',l')] 

+^[(5«(0,2,4;0",1'0-5^(0,2,4;0^2'0-S0(0,2,4;1^2'0 

-5^(0, 2, 4; 0', 3') + 5^(0, 2, 4; O', 2') + 5^(0, 2, 4; 2', 3')) o £3^(0', l', 3') 
+ (5^(1, 2,3;0", 1") -5^(1, 2,3;0",2") -5^(1,2,3; 1",2") 
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-5^(1,2,3; l',2') +5^(1,2,3; l',3') +5^(1, 2,3;2',3'))o£3^l',0',2') 
+ (5§(0, 2, 4; 1", 1") - 5S(0, 2,4; O", l") - 5§(0, 2,4; O", l") 
-S^(0,2,4;0',3') + S^(0,2,4;0',l') + 5^(0,2,4;l',3'))o£0(3',2',0') 
+ (50(1,2,3;1",2")-S«(l,2,3;0",n-S^(l,2,3;0",2") 

-5^(1, 2,3; l',2') + 5«(l,2,3;0',l') + 5^(1, 2,3;0',2'))o£30(2',3',l')]}o 



SI 



(0) 



#5- 



In these subtraction terms the soft eikonal factors contain apart from the soft gluon momentum 
only momenta, which are not affected by the phase space mapping. The soft subtraction term 
da^^^' is only singular in the single soft limit, in all other limits it is integrable. For the integrated 
version we need therefore the integral over the soft eikonal factor for the case where the hard 
partons in the eikonal factor are independent of the phase space mapping. This integral is given 
by [38] 



^3 (Pm Pk,Pa,Pb) = J d(^unres,ijk 5? («, j, k; a, b) 



(141) 



7t 



where 



X = 



4r - -ln(x) -Li2 



X 



+ 0(8), 



{'ai+'al) i^bi + ^bk)' 

Integrated over a one-particle unresolved phase space one obtains 



(142) 



/ 



Agl 
4 2 



(143) 



[Z)0(0, 1,2) (5?(0,2;0', l') -5?(0,2;0',2') -5?(0,2; l',2') 

-53^0,2;0,2)+50(0,2;0,3)+4(0,2;2,3)) 

+dO(0,2, 1) (53^(0, 1;0', 10 (0, l;O',20 -530(0, 1; l',2') 

-53^0, 1;0, 1) +53^(0, 1;0,3) +53^(0, 1; 1,3)) 

+dO(3,2,1)(53°(1,3;1',20-53"(1,3;O', 10 -53^(1, 3;0',2') 

-53^(1,3; 1,3) +53^(1, 3;0,1)+53°(1,3;0,3)) 

+^0(3,1,2) (50(2,3; l',20-53«(2,3;O0 10 -53°(2,3;O',20 

-53«(2,3;2,3)+53°(2,3;0,2)+53«(2,3;0,3))] 

-(t^ + ^) [^3(0,1,3) (53>,3;0',20-5«(0,3;0',10-53>,3;1',20 
-53° (0, 3; 0, 3) + 53°(0, 3; 0, 2) + 53^0, 3; 2, 3)) 
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+A^(0,2,3)(53^0,3;0',2')-53°(0,3;0',1')-53°(0>3;1',2') 
-J3«(0,3;0,3)+530(0,3;0,l) +53^(0,3; 1,3))] }o a 



(0) 



J{|)4. 



4 2 



(144) 



1 Nf-l 

- + — 

4 2 



[£30(0,1,2) (53°(0,2;0',l')-53°(0,2;0',2')-53'(0,2;l',2') 



-53O (0, 2; 0, 2) + 53^(0, 2; 0, 3) + 530(0, 2; 2, 3)) 
+£30(1,0,3)(530(1,3;0', 10 -530(1, 3;0',2') -530(1, 3;1',20 
-530(1,3; 1,3) +530(1,3; 1,2) +530(1, 3;2,3)) 
+£30(2, 3, 0) (530(0, 2; 1', 2') - 530(0, 2; O', l') - 530(0, 2; O', 2') 
-530(0,2;0,2)+530(0,2;0,l) +530(0,2; 1,2)) 
+£30(3,2,1) (530(l,3;l',2')-530(l,3;0',lO-530(l,3;0',2') 
-530(1,3; 1,3) +530(1,3;0,1) +530(1, 3;0,3))] 

+i [£30(0,1,3) (5?(0,3;0',l')-5?(0,3;0',2')-5?(0,3;l',2') 

-530(0,3;0,3)+530(0,3;0,2)+530(0,3;2,3)) 
+£30(1,0,2) (530(1,2;OMO-530(1,2;0',2')-530(1,2;1',20 
-530(1,2;1,2) +530(1,2; 1,3) +530(1, 2;2,3)) 
+£30(3,2,0) (530(0,3; l',2')-530(0,3;0',lO-530(0,3;0',20 
-530(0,3;0,3)+530(0,3;0,l) +530(0,3; 1,3)) 
+£30(2,3, l)(530(l,2;l',2')-530(l,2;0',lO -530(1, 2;0',2') 

-5?(l,2;l,2) + 5?(l,2;0,l)+53°(l,2;0,2))]}o 



(0) 



<i(|)4. 



These integrated contributions cancel all explicit poles of the form as in eq. (130) point by point 
in the four parton phase space. Only in the colour factor NfNc remain some explicit poles, which 
are of the form as in eq. (128). The contribution in eq. (128) vanishes to all order in 8 after the 
integration over the unresolved phase space in D dimensions. 



5 Infrared finiteness 

In this section I show that the three-parton, four-parton and five-parton contributions to the 
NNLO correction are individually infrared finite. There are two types of potential divergences 
which need to be checked: 

(i) Explicit poles in the dimensional regularisation parameter £ = (4 — D)/2, which occur in in- 
dividual stages after one or two integrations have been performed analytically. 
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(ii) Singularities of the phase-space integration. 



We begin with the three-parton contribution. For this contribution only the cancellation of the 
explicit poles in £ needs to be checked. It is easily verified that the combination 

I of o (jaf^+Jaf«^'^ + Jaf ^ + Jaf"'^"'^^) (145) 

is free of explicit poles in £. 



For the four-parton contribution we have to check the explicit poles and the phase-space sin- 
gularities. The four-parton contribution is given by 

I [of [do['^ + daf^''^ - oP o (da'^^"P + rfaf "^"^'^ - daf - daf + Jaf '•"^^'^) 



The combination 



I of (jof +Jar^"^') (146) 

is free of explicit poles in £. This is just the cancellation of the explicit poles in £ of NLO 
calculations. On the other hand as discussed in detail in section 4.7 the explicit poles of the 
combination 

I o (da^^'P + daP'"'''*'' - daf'"''' - daf + Jaf ^^«'^^) (147) 

cancel point by point of the four-particle phase space for the colour factors A^^, A^^, A^^^, Nf/Nc 
and Nj. In the colour factor NfNc the explicit poles in £ cancel after integration over a one- 
particle unresolved phase space in D dimensions. This is related to eq. (128). Alternatively 
they cancel after symmetrisation in the quark momenta. Due to this anti-symmetry finite terms 
obtained from expressing the D-dimensional unresolved phase space measure as an effective 
four-dimensional unresolved phase space measure integrate to zero and can be dropped in the 
numerical program. It is possible to obtain also for the colour factor NfNc a cancellation of the 
poles point by point in the four-particle phase space at the expense of introducing additional 
subtraction terms. However a cancellation of the poles after integration of the one-particle unre- 
solved phase space is sufficient. 

We now consider the phase-space singularities of the four-parton contribution. The following 
combinations are individually integrable over the D-dimensional four-parton phase space: 

of Jaf -of oJa^''^), 

I (of Jaf - oP odaf'-""'^^ , 

I oPo(daP'°'^"''-da^'^''), 

j o^odaf. (148) 
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We now consider the five-parton contribution. Here only phase-space singularities have to be 
considered. The five-parton contribution is given by 



/ 



dof^ - of ) o daf'^' -0^0 (daf '''''' + daf"""'' + daf - daf'"'"^^ ] . (149) 



The potential double-unresolved singularities are: double soft, soft-coUinear, triple coUinear and 
two pairs of coUinear particles. The single-unresolved cases, which need to be checked are one 

soft particle and one pair of coUinear particles. The five-parton contribution consists of two 
channels, the <7gg^^-channel and the qq' gq' q-channel. It is sufficient to check one representative 
case for each type of singularity. We have the following cases: 

• Double unresolved, 0^1^2g3g4^-channel: 

- Double soft: 1^, 2g soft. 

- Soft-coUinear: 1^ soft, 2^ ||3g coUinear. 

- Soft-coUinear: Ig soft, 3^1 1 4^ coUinear. 

- Triple coUinear: 0^| | lg| |2g. 

- Triple coUinear: lg||2g||3g. 

- Two coUinear pairs: 0^| 1 1^ and 2^1 |3g. 

- Two coUinear pairs: 0^| | Ig and 3g| |4^. 

• Single unresolved, 0^1g2g3g4^-channel: 

- Soft: Ig soft. 

- CoUinear: Og||lg. 

- CoUinear: lg||2g. 

• Double unresolved, 0^1^2g3^'4^-channel: 

- Double soft: 1^/, 3qi soft. 

- Soft-coUinear: 2g soft, 1^'||3^' coUinear. 

- Triple coUinear: 1 ^/ 1 1 2g 1 1 3^/ . 

- Triple coUinear: 0^ 1 1 1 1 1 3^' . 

- Two coUinear pairs: 0^| |2g and 1^'| |3^'. 

• Single unresolved, 0^1^'2g3^4^-channel: 

- Soft: 2g soft. 

- CoUinear: 1^'||3^'. 

- CoUinear: l/;'||2p. 
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In the double unresolved limits we have 

lim oi^^ rfaf - oi^^ o dat'*'''' - of ^ o daf"""" = integrable, 

double unresolved -> -> J J 

lim o daf^'' - oP o daf = integrable, 

double unresolved 

lim of^oJa™-^' = integrable. (150) 

double unresolved 

For the single unresolved it is convenient to split the NLO subtraction term da^^"^^^^ into 

^^ii'ng/e ^^ingle, relevant _^ ^^single, remainder (151) 

where do^^^'^sie.relevant ^Qj^^^jj^g jj^g antenna subtraction terms singular in the single unresolved 
limit under consideration and d^^^^^sle, remainder ^Qj^^^jj^^ ^j^g remaining terms. Although the an- 
tenna terms of ^Q^^^"8le,remainder flnite in this particular limit, the matrix element multiplying 
these antenna functions can become singular in an NNLO calculation. (In an NLO calculation 
the infrared-safe observable ensures that there is no contribution from these regions to the NLO 
calculation.) In a similar way we split da'^^^'"^'^ 

^^iterated ^^iterated,relevant _^ ^^iterated,remainder (152) 

such that da!^"-""'^'''^'"'"' is an approximation to daf'^'^''^'"''"' and dal^^rated, remainder ^ 
proximation to dQi^^^^sle, remainder ^ ^asiyQ in single unresolved limits 

lim of Jof - of o da^^sieMevant ^ integrable, (153) 

single unresolved 

lim of ^ O dat'le,remainder _ ^ (3) ^ ^Uerated, remainder ^ ij^^gg^^^le, 
single unresolved 



^ of o (^^ccf - t/af "^'^ - dai^^"^'' - daf*) = integrable. 
The explicit results for the non-trivial singular limits are listed in the appendix B. 



lim 

single unresolved 



6 Conclusions 

In this paper I gave detailed information on the structure of the infrared singularities for the pro- 
cess e'^e~ 3 jets at next-to-next-to-leading order in perturbation theory. Particular emphasis 
was put on singularities associated to soft gluons. The knowledge how to disentangle singular- 
ities at NNLO for this process will also be useful for other processes with three or more hard 
coloured partons, if calculated at NNLO. In this paper a complete list of the subtraction terms 
used in the numerical program "Mercutio2" was given. 
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A Alternative soft subtraction term 



In this appendix I present an alternative form for the soft subtraction terms. The starting point is 
eq. (130) 



/ = X^{a,i,b) \s^{ss!,)-S^isat)-S^{saj) + S^{saj)-S^{si,j)+S^{sbj) 



(0) 



#4, 

(154) 



which in the frame defined by eq. (133) 

1 



Pa 

Pb 

Pa 

Pb 
Pj 



/^(l, 0,0,1), 
/^(l, 0,0,-1), 



2 
1 

0,0,1), 

( 1 , sin 0^ cos (|) , sin 6ft sin (|) , cos 0^) , 
Ej{l, sin 0; , 0, cos Qj) , 



(155) 



leaves uncancelled explicit poles in 1/e for Qj > 0^. The condition Qj > 0^ is equivalent to 
cos Qj < cos 0ft or, expressed in invariants. 



2paPj IpaPb 
> 



^PaibP j '^PaibPb 



(156) 



Let us define two new momenta pa and pj, which are related to and p^ by the Catani- 

Seymour mapping with particle a being the spectator and particle b being the emitter: 



Pa 



^aib 



^ai ~l~ ^ab 

As subtraction term to eq. (154) we can use 



Pa, Pb = Pb + Pi- 



Sib 



Pa- 



A = X^{a,i,b)@{Bj>Bb)[s^{s,y)-S^{sab)-S^{sj},)+S^{saj) 
We then have 



(0) 



I-A = ^X^{a,i,b)- 
(0; < 0ft) In 



(157) 



#4- (158) 



(159) 



' 1 + cos j) ( 1 — cos Oft) 



+0 {Qj > 0ft) In 



2(1 — cos 0ft cos 0y — sin 0ft sin 0j cos (|)) 

(1 -cos0j)(l + cos Oft) 
2(1— cos 0ft cos Qj — sin 0ft sin Qj cos ([)) 



(0) 



<i(|)4 + o(e). 



Integration over the azimuthal angle now leads to zero for the pole terms proportional to 1 /e in 
both cases 0^ < 0ft and 0^ > 0ft. We then have to lift eq. (158) to the five-parton phase space. 
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The terms with S^isab) and S^isaj) pose no problem. The terms with ^3 (5^^) and S^isj},) are 
more subtle. Let us first consider the term 53(5^^). Lifting 53 gives an eikonal factor which 
depends apart from the momentum of the soft gluon on two hard momenta. We would like these 
two hard momenta to be such that under the momentum mapping they are mapped onto pj and 
pg. We can do this by a momentum mapping, which maps five momenta to four momenta. Let us 
denote the five momenta by puPj,Pk,PhPm^ with pj being the soft particle. The eikonal factor 
is given by 

S%j,n) = Sna,^, (160) 

SijSjn 

where pn is a linear combination of pk, pi and pm'. 

Pn = Pk + Pl--. Pm, y2 = • (161) 

i^-yi Skim 

The 5 — > 4-mapping is given by 

Pi = Pi + Pj - z Pn, Pk = Pk, 

l-yi l-yi 

1 i —yi —yi 

Pi = -j Ph Pm = jz rr- rPk, (162) 

l-yi (l-yi)(l-3'2) 



with 



yi = (163) 

Sijm 



Pn is transformed to 



Pn = -r~ Pn= Pk+Pl- T^^Pm, 3^2 = — = 7^^- (164) 

Integration over the unresolved phase space yields 

S-^l?^ j d^unresSl{iJ,n) = S^is-f^), (165) 

with 53 (5) given by eq. (131). 

To lift the term 53(5'^^) we consider four momenta Pi,Pj,Pk,Pi and map them to three momenta 
Pj, Pk and pi. In this case we denote by pi the soft particle. The 4 3-mapping is given by 

Pj = Pi + Pj-yiPn, 

a 

Pk = Pk + yiPn-z — Pi, 
I— a 

Pi = -r—pi, (166) 

I— a 
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with 



Note that 



Pn=Pi+Pj+Pk-7, Ph yi = • a= ■ ■ . (167) 

■^PlPijk mSjn yiSkn+yiSln + Skl 



With 



the eikonal factor is given by 



Pn = Pj + Pk-7^ Pi- (168) 

^PiPfk 



Pn = (l-yi)Pn, (169) 



5^(7, = 87ia,^. (170) 

S jiSin 



Integration over the unresolved phase space yields 

S-^ll^ j d(Sfunres^^{jM) = S^{sjf,). (171) 

We remark that the subtraction term of eq. (158) treats the partons a and b asymmetric for the 
construction of the hatted momenta, with a being the spectator and b being the emitter. A more 
symmetric version is given by 

K^'^aii'^ib y-^FaihFj ^Faibfb/ '- 



-^3 



(0) 2 



Finally we remark that the expressions of eq. (154) and eq. (158) give rise to finite terms, which 
have to be taken into account. The first source of these terms is straightforward and arises from 
the expansion of the functions 53 to order which occur in eq. (154) and eq. (158). The 
terms of the function Z3 (a, /, b) give no contribution to the finite terms, they vanish for the same 
reasons as the pole terms. However the D-dimensional phase space measure d^s^unres gives an 
additional contribution at order £*^. Let us consider again the frame defined in eq. (133) and in 
eq. (155). In four dimensions we parametierised by 

Pb = £'i(l,sin0icos(|),sin0^sin(|),cos0i). (172) 

Qt) is the polar angle and ^ the azimuthal angle. In Z)-dimensions we can parameterise pt, by 

Ph = £'/,(l,sineicos02,sineisine2COS03,cos0i,...). (173) 

The dots indicate the components in the additional dimensions. 0i, 02 and 63 are polar angles. 
As in the chosen frame all other particles lie in the jc — z plane, nothing depends on 03 and we 
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can integrate over 63 and all additional angles 64, Qd-2- The unresolved phase space in 
D-dimensions is then given by 



J d^unres,aib = {xiX2X3Saib) ^ J d^X d{l - Xi - X2 - X3) 

[ Je2(sine2)-2^ (174) 



where 



^ai ^ib ^ab 
X\ = , X2 = , X3 = . (175) 

^aib ^aib ^aib 

We would like to write this as an effective four-dimensional phase space. For an integrand which 
depends on the polar angle 62 only through cos 62 we may write 

2jdQ2isinQ2y^ficosQ2) = J d^\sm^\'^f(cos^). (176) 


This allows us to treat the polar angle 62 of the D-dimensional parameterisation as the az- 
imuthal angle (|) of the four-dimensional parameterisation. We can therefore write the unresolved 
phase space in D-dimensions effectively as a phase space in four dimensions with additional 
8-dependent factors: 

J d(^unres,aib = (47t)^ ixiX2X3S aib)'^ J d^X d{l - Xi - X2 - X3) 



2ji 



d^\sm(^\-^. (Ill) 



Expansion of | sin(|)|^^^ to order £^ gives a finite contribution when combined with the l/£ poles 
of 53. The e^-terms of the factors {'^'Ji)^{xiX2X3Saib)~^ and r(l/2)/r(l/2 — e) can be dropped, 
as they vanish after integration over the unresolved phase space for the same reason as the pole 
terms. 



B Singular limits of the five parton contribution 

In this appendix limit 1 list the explicit forms of the non-trivial singular limits of the five parton 
contributions. The double unresolved limits of do^^^ have been considered in [27,93-98], the 
singular limits of the four-parton tree level antenna functions have been taken from [37]. The 
singular functions are denoted as follows: The single soft factor is denoted by Sijk, where j is 
a soft particle between particles i and k. The double soft factor is denoted by Sijki, in this case 
particles j and k are soft. The singular function for two coUinear particles i and j is denoted by 
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Pij. No distinction is made in the notation between g — > gg-, g — > qq- and q — > ^g-splittings. This 
can be deduced from the involved particles. If particles / and j become coUinear the resulting 
momentum is denoted by {ij). The singular function for the triple coUinear limit is denoted by 
Pijk- Finally the soft-coUinear function is denoted by SCiji^it^i-^ with j the soft particle and k and / 
the coUinear pair. We note that in the soft-coUinear limit the following combinations 

{SCij(^ki) - Sijk) Pkh 
{SCij{ki) - SCiji^ik)) Pkh 

{Sim-Sijk)Pkl (178) 

are all integrable. 

In the double unresolved limit we have for the 0^1g2g3g4^-channel: 



Double soft: Ig, 2g soft 



lim da^^^ = 7 I ~r ["^0123 + 5'o2i3 + 5'oi35'324 + S023S314 + 53124 + 53214] 
1,2^0 ^ 6 [ 4 

+ ^ [-5'oi24 - 5'0214 — 5'0135'024 — 5'3145'o24 - 5'oi45'o23 " 5'oi45'324 + 5'oi45'o24] 



1 



4N^ 



50145024 ? o 



(0) 



#5, 



/ \ 1 r 

1^2% ("^^3"" ^) ~ 6\~4^ ["'^0123 - 5'0213 - 5'3124 " 53214 " 5oi35o23 " 53145324 + 5oi45o24] 



+ ^ [5oi24 + 5o214 +5oi45o24] + [— 5oi45o24] } ° 



(0) 



#5, 



1^2% ~ 6 { [—"^0135324 — 5o235314 + 5oi35o23 + 53145324 — 5oi45o24] 



(0) 



#5- 



+ ^ [5'0135'024 + 5'0235'014 + 5'3145024 + 53245014 -25oi45o24] J o 

Soft-coUinear: 1^ soft, 2g\\3g coUinear. 

^Jim Jaf = [^Coi(23)/'23+5Coi(32)/'23 + 5C4i(32)/'23 + 5C4i(23)/'23] 



(179) 



+ ^[-25014^23] jo ^f^'^ 



#5, 



^Jm^^^ (-rfaf "^'^) = ^ 1^ [-5Coi(23)/'23 -5Coi(32)^23 -5C4i(23)/'23 -5C4i(32)^23] } 



(0) 



#5, 



Jim^^(-<,af-)=i{il25„,,P.3]}o 



(0) 



(180) 
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Soft-coUinear: Ig soft, 3p\\4o coUinear. 



M _ 1 / 



jjim^^^rfa^ = 6 [f ['5012^34 + 5C2i(34)P34] 

+ ^ [-'^012^34 -'5Coi(34)'f34-'SC2i(43)/'34] + ^^'^Ql(43)'P34 



(0) 



#5, 



lim 

1-^0, 3||4 



-da: 



double 



\ 1 / 
J ( 



[-5C21 (34)^34 - 5C21 (43)^34 + SQi (43)^34] 



"4 ['^Ql(34)'P34+5'Coi(43)i'34] + [-'5Ql(43)'f34] 



(0) 



lJi)"^||4("'^"'"""0 = [-5oi2P34 + S2i(34)^34-Soi(34)/'34] 



+ T [^'012^34 + 5'21(34)^34 - ^'01(34)^34] 



(0) 



#5- 



Triple coUinear: 0^| 1 1^1 |2g. 



o|f}|f/^f = ^ ^ [^012+^021] + 7 [-^^012-^021 -P^'n] + 7^ [^0 
lim 

0||1||2 



1 (N} 



1 



1 



4N^ 



012J 



(0) 



(181) 



lim (-Jaf«^'^) = i|^[-P012-P02l] + ^[P012 + P021+P^f2] 



4A^ 



012] 



(0) 



#5, 



lim (-daf"""") = integrable. 

0IIIII2V J 

Triple coUinear: lg||2g||3g. 



(182) 



|m ^/af = ^|^[Pl23+ ^132 + ^213 +/'231+ ^312+^321] jo 

lim (-daf'*''^') = i-Pm - Pm - Pin - P231 - P312 - P321] ] 

\\2\\3\ ' J 6{ 4 ' 



1||2||3 

lim (-da"''""''^ 

1\\2\\3\ 



(0) 



integrable. 



#5, 
(183) 



Two coUinear pairs: 0^| 1 1^ and 2g| |3g. 



lim do 

0||1,2||3 ■ 

lim 

0||1,2||3 

li 

0||1 



(0) 



^ I ^ [P01P23 +P01P32] + ^ [-P01P23 - P0lP32] 



}• 



(0) 



#5, 



lim i-daf 



^[-^01^23-^01^32] 



(0) 



S|3 (-^<"") = l{\[P0lP23+P0lP32]]o 



SI 



(0) 



#5, 



#5- 



(184) 
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Two coUinear pairs: 0^ 1 1 Ig and 3g\ \4q. 



lim /of - ^ l^[P0lP34] + \ [-P01P43-P10P34] + t4o A)1^34 



0||1,3||4 

lim 

0||1,3||4 



1 (N} 



1 



4N} 
1 



(0) 



= Z i ^ [^0l/'34] + 7 [2^01^34] + 4^ [-P01P34] 



(0) 



lim f-Jaf'"''^'') 

111.3114 V / 



1^ [2P01P34] !> o 



(0) 



#5- 



0||1,3||4V - / 6 

In the single unresolved limit we have for the 0^1g2g3^4^-channel: 



Soft: Ig soft. 



#5, 
(185) 



lim ^-daf""^) = ^ { ^ [" (^3(0' 2, 3) +D«(0, 3, 2)) (5oi2 + Sou + S213) 

-(D^(4,3,2)+D^(4,2,3)) (^214 + ^314 + ^213) + (A^(0,2,4) +A^(0,3,4)) ^014] 
+^[a0(0,2,4) (%2 + 52i4)+A^(0,3,4) (%3 + ^314) + (A^ (0,2,4) +aO(0, 3,4)) Sqm] 



[- (a0(0,2,4) +a0(0,3,4)) %4] | o 



(0) 



#5, 



lim dC^^ratedMevant ^ 1 | ^ ^^^O^q^ 2, 3) +^^(4, 3, 2)) (S012 + S213 + S314) 

+ (D§ (0, 3, 2) + D§(4, 2, 3)) {Sou + S213 + S214)] 

+^ [- (A«(0,2,4)+a0(0,3,4)) (5oi2 + 5013 + ^214 + 5314) 

-(d0(0,2,3)+d0(0,3,2)+d0(4,3,2)+d0(4,2,3)-A0(0,2,4)-A0(0,3,4))%4] 



4N^ 



[(a0(0,2,4)+a0(0,3,4))%4] \o 



SI 



(0) 



#5, 



lim (- Ja?--) = 1 { ^ 1 [D«(0, 2, 3) (^013 - ^314 + ^OM + 5^132 - %14 - ^^14) 

+Z)0(0, 3, 2) (5oi2 - 5214 + 5oi4 + ^0 
+^3(4, 3, 2) (S214 — 5'oi2 + 'S'oi4+5'^^^ 



'03123 '^2314 "^0314) 



+D^(4,2,3) (5314-5013 +%4+ 5: 



32142 



"^0123 0143 
"^0132 "^0142 



) 



+A3 (0, 2, 4) (-5oi3 - 5314 - 5oi4 - 5o2i^ + 5q213 + 55413) 

+^3(0, 3,4) (-5oi2 - 5214 - 5oi4 - 5q3J43 + 5q3J2 + 534J2)] 

+^[(dJ(0,2,3)+Z)§(0,3,2)+D?(4,2,3)+D§(4,3,2))5oi4 
+A^(0, 2, 4) (-5oi4 - 552142 + 5q213 + 5^413) 
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(0) 



+A^(0,3,4) (-^014-5o3j43+5^j2 + '^3412)]} 

lim [-daf') = 1 1 ^ 1 [dO(0, 2, 3) (5oi3 - ^314 - ^om - 555^35 + Sj^^^^ + S^^J 

+^3(0, 3, 2) (S012 -5214 -^014 -'S'q3^23 +"^2314 +"^0314) 

4-^3 (4, 3, 2) (5214 - 5'oi2 - 5*014 - 5^^^ + 5^^^ + 5q^^) 
+^3(4, 2, 3) (5314 - 5oi3 - 5oi4 - 5^21^ + 5q^52 + "^0142) 



) 



02142 "^0213 "^2413 
03143 ~ "^^12 ~ "^3412)] 



(0, 2, 4) (5oi3 + 5314 - 5oi4 + 5, 
+A§ (0,3,4) (5oi2 + 5214 - 5014 + 5, 

+^ [aO(0,2,4)(Soi3 + 5314-5oi4 + 5, 

+A3(0,3,4) (5o12 + 5214 -5oi4 + 5^j^-5^j2~ '^3412." ^° 



02142 '^0213 '^2413^ 



(0) 



#5- 



(186) 



CoUinear: 0^||1^. 



lim(-Jaf"^^^) = i|^[(-2D^(01,2,3)-2DO(01,3,2)+AO(01,2,4)+A^(01,3,4))Poi] 



0||1 V 



+^[2(A§(01,2,4)+A§(01,3,4))Poi] 
+iP [- (^3(01,2,4) +A«(01,3,4))Poi] I o 



(0) 



#5, 



^ 6 



1 {N} 



0||1 



[(dO(01,2,3)+dO(01,3,2)+dO(4,3,2)+D«(4,2,3))Poi] 



1 



+- [-(d0(01,2,3)+d0(01,3,2)+d0(4,3,2)+d0(4,2,3)+A«(01,2,4) 



(0) 



#5, 



+A§(01,3,4))Poi]+^[(A§(01,2,4)+A§(01,3,4))Poi] 
Urn (-Ja?'--^) = [(-D0(4,3,2)-D§(4,2,3)+d0(01,2,3)+d0(01,3,2) 
-AO(01,2,4)-AO(01,3,4))Poi]+^[(£>3(01>2,3)+dO(01,3,2)+dO(4,2,3) 

+Z)0(4,3,2)-AO(01,2,4)-A^(01,3,4))Poi]} 

lim ( —daf'!'-^^ ] = integrable. 

0111 V / 



(0) 



0||1 

CoUinear: lp||2 



(187) 



1 (N} 



lim [-dat-'^^) = - 1^ [2(-D0(0, 12,3) -D^(0,3, 12) -^0(4, 12,3) -7)0(4,3, 12))Pi2] 
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+^[2A?(0,12,4)Pi2] 



(0) 



lim d^!^^ratedMes>ant (^O^q^ 3) +D^(0, 3, 12) +^0(4, 3, 12) +^0(4, 12,3)) Pxt\ 



1||2 



1 



+^ [-2(A§(0,12,4)+A^(0,3,4))Pi2] 



0/' 



(0) 



#5, 



lim (-Ja^--) = i I i [2A0(0, 3,4)^1^] | o 
lim ( —da^'f'^' ] = integrable. 

1||2 \ / 



(0) 



#5, 



(188) 



In the double unresolved limit we have for the 0^1^/2^3^4^-channel: 



Soft-coUinear: 2g soft, 1^'||3^ coUinear. 



lim da 



(0) 



2^0, 1||3 



T + -A: I I -J- (5'Co2(31)^13 +5'C42(13)^13) 



1 



4Nr 



(-5024^13 - 25Co2(31)^13 + 25Co2(l3)^13 " '2-SC 42(^13) Pl3 + 25C42(3i)Pl3) 



(0) 



#5, 

lim f-rfaf"^'^) 

-^0,1||3V J 



lim 

2^0,l||3 



-da' 



almost^ 



I Nf-l 

4 2 

1 A^f-1 

- + — 

4 2 



(-SCo2(31)^13 -5'C42(i3)i'l3) 



(0) 



#5, 



A^c 1 

-^-^ (5'021^13 -5'023^13 + 5'423^13 -^'421^13) 



1 



4A^c 

Triple coUinear: 1^/ 1 |2g| |3^/ 



(5024^13 - 25021^13 + 2^023^13 - 2^324^13 + '^SmPu) \ ° 



SI 



(0) 



#5- (189) 



lim d<3 

1||2||3 ■ 



(0) 



1 I Nf-\ 
4 2 



^(P213+^23l) + 4^(-Al3)^ 



SI 



(0) 



#5, 



lim (-daf-''^\ = T- + 

iii2||3v ^ y \^4 2 

lim i-da^^'^A = 0. 

iii2||3V y 

Triple coUinear: 0^| 1 1^ 1 13^/. 



Nr 



1 



123 



(0) 



#5, 



(190) 



lim da 

0||1||3 



(0) 



+ ^^1^013 + 77^(7 + ^:^ ) (-^013) 



4 V4 



4Nr V4 
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id 
013 



+ 



1 fl 



AN} V 2 



id 
013 



SI 



(0) 



lim (-(^a' 

0||1||3V 

1 /I 

+ 



double \ 



Ncfl Nf-l 
4 U 2 



(-/'013) + 



1 /I % - 1 



4A^r V4 



- + 



^013 



4 V2 



013 



+ 



1 



4N} 



1 



old 
013 



(0) 



lim (-da^'"''"') = 0. 

0||1||3V J 
Two coUinear pairs: 0^| |2g and 1^'| |3^'. 



(191) 



lim do 



(0) 



0||2,1||3 

lim f-Jaf"^'^) 

0||2,1||3V J 

lim f-Jaf'«^'0 

112. 1113 V / 



1 Nf-\ 

- + — 

4 2 



^(^^02^13) + ^ (-^02^13) 



(0) 



Nc (I Nf-l 



0||2,1||3 



4 V4 

1 Nf-l 

- + — 

4 2 



(-^02^13) ^ o 



(0) 



4A^. 



(^02^13) ^ o 



(0) 



#5, 



#5- 



(192) 



In the single unresolved limit we have for the 0^1^'2g3^'4^-channel: 



Soft: 2g soft. 



4 



1 I 

4 2 



(-2(^0(0,1, 3)+£30(3,4,0))5o23 



-2(£?(l,0,4)+£?(4,3,l))5i24) 

+^(-2(£30(0,l,4)+£0(4,3,0))5o24-2(£3«(l,0,3)+£30(3,4,l))5i23) 



4Nr 



I Nf-l 



((£0(0, 1,3) +£30(4,3, 1)) 5i23 + (E^ihOA) +£3^(3,4, 0)) 5o24) 



+i ((£0(0, 1,4) +£0(3,4, l))Si24 + (£3'(l:0:3)+£30(4,3,0))5o23) 



(0) 



d(^5, 



iterated 



Nr 



^ + ^^^ ) (£3(0:l:3)+£30(4,3,l)+£30(l,0,4)+£30(3,4,0))(5i24 + So23) 



+ ^(£30(1,0,3) +£30(4,3,0) +£30(0, 1,4) +£30(3,4, l))(5o24 + 5i23) 



4iV. 



1 I Nf-l 
4 2 



(-£30(0, 1,3) -£30(4,3,1) -£30(1,0,4) -£30(3,4,0)) 



(25o23 + 25i24 — 25o21 — 25324 + ^024 +5123) 

+^(-£30(1,0,3) -£30(4,3,0) -£30(0, 1,4) -£30(3,4,1)) 
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(25o24 + 25i23 - 25021 - 2^324 + 5'o23 +-$'124)]} ° 



(0) 



lim 

2^0 



3 ,qci' gq' q-^qggq-*qgq " ""3 ,qq' gq' q-*q(f c/ Q-^QgQ 



^almost 



)- 



4 2 



1 , Nf-\ 
4 2 



[£30(0, 1 , 3) (5o23 - 2Si24 + S024 + ^324 + 5, 



01231 "^4231 '^0124) 



+£■3(1, 0,4) (5i24 - 25o23 + 5'l23 + 5'324 + "^10240 ~ "^3240 ~ "^1023 
+£3 (4, 3, 1) (5i24 - 25o23 + 5'024 + ^021 + 5'^2T3 ~ "^0213 ~ "^02^) 

34204 ~ "^1204 ~ "^1234 



) 



+£■3 (3,4, 0) (5o23 - 25i24 + 'S'123 + 5'021 + S:^2QA ~ "^1204 ~ "^1234)] 

+ i [£30(0,1,4) (5o24-25i23+So23+5324 + 5, 

+£3 ( 1 , 0, 3) (5123 - 2^024 + 5124 + ^324 + 'S'io230 ~ "^4230 ~ "^1024) 
+£?(3,4, 1) (5i23 -25o24 + 5o23 +5021 +5: 



01241 "^3241 "^0123^ 



34214 "^0214 "^0234) 
+£3 (4, 3, 0) (^024 - 2^123 + 5124 + 5o21 + ^^2^ ~ "^1203 ~ "^1243)] ] 



AN, 



1 Nf- 
4 2 



[(£0(0, 1,3) +£3^(4,3, l))5o24+ (£3"(1:0:4) +£0(3,4,0)) 5i23 



+2(£30(0, 1,3) +£30(4,3, 1) + £30(1, 0,4) +£30(3,4, 0))(5o23+5i24-5o2i -5324)] 
+^ [(£30(0, 1,4) +£30(3,4, 1)) 5o23 + (£3 (1,0,3) +£30(4,3,0)) 5i24 

+2 (£30(1,0,3) +£30(4,3,0)+£30(0, 1,4) +£30(3,4, l))(So24 + 5i23 -5021-5324)]]} 



(0) 



#5, 



A^l 
4 2 



1 Nf- 
4 2 



[£3(0, 1,3) (5o23 — 5o24 — 5324 — 5| 



&0 ("^''5^^.' J 

+54231 +5qj24) +^3(1>0>4) (5l24 -5i23 -5324 -5^q24o+'53240 + "^1023) 
+£3 (4, 3 , 1 ) (5i24 - 5o24 - 5o21 - 543213 + 5q213 + ^'0243) 
+£3(3,4,0) (5o23-5l23-5o21 -"^34204 +'^1264 +'^1234)] 
+ 4 [^3 (0, 1 , 4) (5024 - 5o23 - 5324 " 5^241 + "^3241 + "^0123) 
+£3(1 , 0, 3) (5i23 - 5i24 - 5324 - 5^0230 + "^4230 + "^1024) 
+£3(3,4,1) (5i23-5o23-5o21 -'^34214 + "^0214 + "^0234) 
+£3(4,3,0) (5o24-5i24-5o21 ~ -^43203 +'^1203 +"^1243)] } 

CoUinear: 1^/||3^'. 
In this case we have 

Jterated 



01231 



(0) 



(193) 



^^4,qq'gg'q^qg'q'q-^^3,qq'gq'q^qq'q'q^qgq " mtegrable. 



'^^Ifq'gq'q = intcgrablc. 



(194) 
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The non-trivial limits are: 



Um(-^aS'.) = Q + (-^^(0,13,2)-D0(0,2,13)-D0(4,13,2) 



-D0(4, 2, 13)) Pn + 4^^3(0, 13,4)Pi3 1> o 



SI 



(0) 



#5, 



\_ Nf-\ 

4 2 



(Z)0(0, 13,2) +d0(0,2, 13) +^0(4, 13,2) 



+D«(4, 2, 13)) + ^ (-A0(0, 13, 4) - A«(0, 2, 4)) Pi3 



(0) 



#5, 



'™(-''<X'. 



(0) 



qq' gq' q-^qggq^qgq 
2 



j^almost 

'iM'gq'q^qq'q'q-'qgq 



#5- 



(195) 



CoUinear: l^||2g. 



dct 



single 

''^M' gq' q^qggq 



-da 



iterated 

"iMg^q-^qggq^qgq 

^^3,qq'gq'q 



integrable, 
integrable. 
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lim ( -da' 

1||2 



double 

^:qq'g<fq^qq'q'q-*qgq 



^ +^^^ ) (-2£30(12,0,4)-2£30(4,3,12))Pi2 



+^ (-2£30(12, 0, 3) - 2E?(3, 4, 12)) P12 



+ 



1 



4Nr 



I Nf-l 
4 2 



+£30(4,3, 12))Pi2 + - {E^{0, 12,4) +£3^(3,4, 12))Pi2 



(0) 



(^3(0,12,3) 



#5, 



lim ^a^terated 

nm aa3.afl'pfl' 



1112 



3,qq'gq'q^qq'q'q^qgq 



Nr 



1 A^f-1 



4 + ' 2 



(£?(0,12,3)+E?(4,3,12) + £30(12,0,4)+£?(3,4,0))Pi2 



+^(£30(12,0,3) + £:3V,3,0) + £30(0,12,4)+£30(3,4,12))Pi2 



4Nr 



1 A^f-1 



4 + - 2 



(-£0(0, 12,3) -£3° (4, 3, 12) -£0(12,0,4) -£0(3,4,0)) P12 



+^ (-£30(12,0,3) -£30(4,3,0) -£30(0, 12,4) -£30(3,4, 12))Pi2 



(0) 



#5, 



lim ( - Ja' 

1112 V 



1||2 



almost 

'^A^gq'q^qggq^qgq 



dal 



almost 



4 



1 A^f-1 

- + — 

4 2 



(£30(12,0,4) -£30(3,4,0) -£30(0, 12,3) +£30(4,3,12))Pi2 
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+\ (£0(12,0,3) -£0(4,3,0) -£0(0, 12,4) +£0(3,4, 12))Pi2 



4 
1 

^4A^ 



i + (^3 (12,0,4) +£0(3,4, 0))Pi2 



+i(£30(12,0,3) + £30(4,3,0))Pi2 



-^3 



(0) 2 



^^^5- (197) 
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